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Abstract 

We study the coupling of N charged scalar particles plus the electro-magnetic field to ADM 
tetrad gravity and its canonical formulation in asymptotically Minkowskian space-times without 
super-translations. To regularize the self-energies both the electric charge and the sign of the 
energy of the particles are Grassmann- valued. The introduction of the non-covariant radiation 
gauge allows to reformulate the theory in terms of transverse electro-magnetic fields and to extract 
the generalization of the Coulomb interaction among the particles in the Riemannian instantaneous 
3-spaces of global non-inertial frames, the only ones allowed by the equivalence principle. 

Then we make the canonical transformation to the York canonical basis, where there is a sep- 
aration between the inertial (gauge) variables and the tidal ones inside the gravitational field and 
a special role of the Eulerian observers associated to the 3+1 splitting of space-time. The Dirac 
Hamiltonian is weakly equal to the weak ADM energy. The Hamilton equations in Schwinger 
time gauges are given explicitly. In the York basis they are naturally divided in four sets: a) the 
contracted Bianchi identities; b) the equations for the inertial gauge variables; c) the equations for 
the tidal ones; d) the equations for matter. 

Finally we give the restriction of the Hamilton equations and of the constraints to the family of 
non-harmonic 3- orthogonal gauges, in which the instantaneous Riemannian 3-spaces have a non- 
fixed trace 3 K of the extrinsic curvature but a diagonal 3-metric. The inertial gauge variable 3 if 
(the general-relativistic remnant of the freedom in the clock synchronization convention) gives rise 
to a negative kinetic term in the weak ADM energy vanishing only in the gauges with 3 K = 0: is 
it relevant for dark energy and back-reaction? 

In the second paper there will be the linearization of the theory in these non-harmonic 3- 
orthogonal gauges to obtain Hamiltonian post-Minkowskian gravity with asymptotic Minkowski 
background, non-flat instantaneous 3-spaces and no post-Newtonian expansion. This will allow to 
explore the inertial effects induced by the York time 3 K in non-flat 3-spaces (they do not exist in 
Newtonian gravity) and to check how much dark matter can be explained as an inertial aspect of 
Einstein's general relativity. 



2 



I. INTRODUCTION 



By re-expressing the 4-metric in terms of tetrads inside the ADM action, an ADM formu- 
lation of tetrad gravity was presented in a series of papers [1, 2, 3, 4] for globally hyperbolic 
space-times asymptotically Minkowskian, parallelizable (so to admit ortho-normal tetrads 
and a spinor structure [5]) and without super-translations (see Refs.fl, 3] for the needed 
boundary conditions on the 4-metric and on the tetrads). This allowed the development 
of canonical tetrad gravity with its fourteen Dirac constraints and the identification of the 
York map as a Shanmugadhasan canonical transformation adapted to ten of the constraints 
[6]. As a consequence we now have a York canonical basis in which the tidal effects of the 
gravitational field (the polarization of the gravitational waves in the linearized theory) are 
separated from the inertial effects, which are described by the gauge variables conjugated to 
Dirac first-class constraints. Even if this separation is only 3-covariant on the instantaneous 
non-flat 3-spaces corresponding to a clock synchronization convention, it allows to give a 
physical interpretation of all the quantities appearing in tetrad gravity and to rewrite the 
ADM Hamilton equations in a new form. 

ADM tetrad gravity is formulated in an arbitrary admissible 3+1 splitting of the globally 
hyperbolic space-time, i.e. in a foliation with instantaneous space-like 3-spaces tending to 
a Minkowski space-like hyper-plane at spatial infinity in a direction independent way: they 
correspond to a clock synchronization convention and each one of them can be used as a 
Cauchy surface for field equations. As shown in Refs.fl, 3] the absence of super-translations 
implies that the SPI group of asymptotic symmetries is reduced to the asymptotic ADM 
Poincare' group and the allowed 3+1 splittings must have the instantaneous 3-spaces tending 
to asymptotic special-relativistic Wigner hyper-planes orthogonal to the ADM 4-momentum 
in a direction-independent way. At spatial infinity there are asymptotic inertial observers, 
carrying a flat tetrad e A (^rj^e^e 1 ^ = ^t/ab), whose spatial axes can be identified with the 
fixed stars of star catalogues. 

Radar 4-coordinates a A = (r; a r ) adapted to the 3+1 splitting and centered on a time-like 
observer a+(r) are used instead of local 4-coordinates x M . If a+ h- > a A is the coordinate trans- 
formation to these adapted 4-coordinates, its inverse a A i— > a+ = z^(t, a r ) defines the embed- 
ding of the instantaneous 3-spaces S T into the space-time. The gradients z A (r, a r ) = 9Z q^a ~ 
are the transition functions for transforming tensors: for the 4-metric we have A gAB{ T , < yT ) — 
z^(r,a r ) z v B (r,(j r ) A g tlu (z(T,(? r )), so that we get y / \det A g llu \ = y /\det I g A ~B\/\det z%\. Then 

we put the decomposition 4 gAB( T ,°~ r ) = 4 E A a \r, a r ) 4 ^( a )(/3) 4 E B 8 \r, <r r ), with 4 E A \r, a r ) 
arbitrary cotetrads, inside the ADM action to obtain ADM tetrad gravity. 

The three independent space-like 4- vectors z^(r,cr) are tangent to S r in the 
point (r, cr r ). Instead the r-gradient of the embedding has the standard de- 
composition z%(t, a r ) = ^(1 + n) P + n r z^j(r,a r ) along the unit normal 

z a (t, a r ) l A {r, a r ) and the tangents to S r defining the lapse (N — 1 + n) and shift (N r = n r ) 
functions. 

Let us remark that if we find a solution A gABiji^ r ) of the ADM Hamilton equations in 
some gauge, then it has an associate preferred 3+1 splitting of space-time, whose instan- 
taneous 3-spaces are dynamically determined (dynamical clock synchronization convention) 
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by the given extrinsic curvature and by the given lapse and shift functions, as noted in 
Refs. [7]. As a consequence, there is a preferred world 4-coordinate system adapted to it. 
Let us take the world-line of the time-like observer as origin of the spatial coordinates, i.e. 
rr M (r) = (x°(r);0). Then the space-like surfaces of constant coordinate time x°(t) = const. 
coincide with the dynamical instantaneous 3-spaces S T with r = const, of the solution. 
Then the preferred embedding is 

x » = z i*{t, a r ) = x»(t) +e?a r = 5% x°(r) + e£ a r . (1.1) 

If we choose the asymptotic flat tetrads e A = 5% S A + 5f S A and x°(t) = x° + e° r = x° + r, 
we get z»(T,a r ) = 5% x° + e% z^ A (r,a r ) = and %„(x = z(r,a r )) = e A e^g AB (r,a r ) 
(e A are the inverse flat cotetrads). Then by means of 4-diffeomorphisms we can write the 
solution in an arbitrary world 4-coordinate system in general not adapted to the dynamical 
3+1 splitting. This gives rise to the 4-geometry containing the given solution. 

In this paper we are going to study ADM tetrad gravity in presence of the following 
matter: N positive-energy charged scalar particles plus the electro-magnetic field. 

The isolated system of N positive-energy charged scalar particles (with Grassmann- valued 
electric charges for a UV regularization of self-energies) with mutual Coulomb interaction 
plus a transverse electro-magnetic field in the radiation gauge have been studied in special 
relativity in Refs. [8, 9], where its inertial rest-frame instant form was developed: in it 
the instantaneous 3-spaces (the Wigner hyper-planes) are orthogonal to the conserved 4- 
momentum of the isolated system. The extension to non-inertial rest frames, in which the 
instantaneous 3-spaces S T tend to Wigner hyper-planes at spatial infinity (asymptotically 
orthogonal to the conserved 4- momentum in the non-inertial rest-frame instant form), was 
done in Ref.[10]. 

The starting point of these developments was the parametrized Minkowski theory for 
isolated systems [11, 12] (positive-energy particles, strings, fields, fluids; see also Refs. [13, 
14, 15, 16, 17, 18] ) admitting a Lagrangian formulation, the precursor of the described 
formulation of ADM tetrad gravity. After an admissible 3+1 splitting of Minkowski space- 
time with the instantaneous 3-spaces £ r described by embedding functions x^ = z^{r,a r ) 
[a A = (t; a r ) are radar 4-coordinates centered on an arbitrary time- like observer] one defines 
a Lagrangian on S r depending on the given matter and on the embedding through the 
induced 4-metric 4 qab{t, o~ r ) = (z A n^ z^)(r, a r ), This Lagrangian is obtained by the matter 
Lagrangian by coupling it to an external gravitational field and by replacing the external 4- 
metric with A g A B[z{r, a r )] (a functional of the embedding) after having redefined the matter 
fields so that they know the clock synchronization convention (for a Klein-Gordon field (f>(x) 
we use 0(t, a r ) = (j>(z(r, o" r ))). Each admissible 3+1 splitting corresponds to a convention 
for clock synchronization and defines a global non-inertial frame centered on the observer 
and the components of the 4-metric A gABij, °~ r ) play the role of the inertial potentials in the 
non-inertial frame. 

As said, the Lagrangian of parametrized Minkowski theories depends on the embeddings 
z^(t, a u ) and on matter variables adapted to the foliation. Due to the invariance of the action 
under frame-preserving diffeomorphisms, the embeddings z^ir^ a u ) are gauge variables. This 
implies the gauge equivalence of the description of the isolated system in any non-inertial 
or inertial frame, namely its independence from the clock synchronization convention and 
from the choice of 3-coordinates on S T . 
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The transition to globally hyperbolic space-times is done in accord with the equivalence 
principle, according to which only a description of tetrad gravity plus matter in global non- 
inertial frames is possible (inertial frames exist only locally near a free falling particle). In 
particular the absence of super-translations identifies the non-inertial rest frames (orthogonal 
to the ADM 4-momentum at spatial infinity) as the only relevant ones. 

The real difference with parametrized Minkowski theories is that the independent vari- 
ables are now the cotetrads 4 E A a \r, a) (or the 4- metric 4 gAB( T ,° t ) i n metric gravity) and 
not the embeddings ^(r, a) (z A are only transition coefficients, while in special relativity 
they are flat cotetrads), which are dynamically determined a posteriori instead of remaining 
gauge variables as it happens in special relativity. 

Following Ref. [11, 13] we will describe N scalar particles of masses rrii, i=l,..,N, with 
3-coordinates ?7j(r) on the instantaneous 3-spaces E T (diffeomorphic to R 3 ) identified by the 
intersection with E r of their world-lines a;f(r) = z^(r, ?y[(r)) in M 4 . Therefore the world- 
lines xf(r) are derived quantities (as shown in Ref. [8] they are covariant non-canonical 
predictive coordinates) and, as shown in Ref. [6], describe particles with a definite sign of 
the energy (due to clock synchronization each particle is described by 3, not 4, position 
variables). As a consequence, the momenta pf(r) are not defined: as shown in Ref. [8] it is 
possible to define them so that the mass shell constraints ep 2 — m 2 c 2 « are satisfied also 
in presence of interactions. However other definitions are possible. 

To regularize the gravitational self-energies we will assume that the signs rji of the energies 
of the particles (a topological number with the two values ±1) are described by two complex 
conjugate Grassmann variables 9i, 9*: rji = 9* 9i, rjf = 0, rji rjj = rjj rji ^ for i ^ j (after a 
formal quantization they describe a 2-level system). In this way we are implementing a % ^ j 
rule like it was done in classical electrodynamics by using Grassmann variables to describe 
the electric charges of charged relativistic scalar particles [14, 15] : Qi = 9^* 9^ Q \ Q\ = 0, 
Qi Qj = Qj Qi 7^ for i 7^ j, to avoid divergencies in the electro-magnetic self-energies. 

The electro-magnetic field is described by the vector potential A a (t, <j u ) = 
z^(t, cr u ) A^(z(t, a u )) and by the field strength F ab (t, cr u ) = d A A b (t, cr u ) - d B A a (t, a u ) = 
z a (t,<j u ) z b (t,<j u ) F^ v (z{t,(j u )). Following Ref. [10], we shall restrict the formulation to an 
electro-magnetic field in the radiation gauge in the instantaneous 3-spaces E T . Even if this 
gauge is not covariant, it allows to extract the non-inertial analogue of the Coulomb potential 
among the particles and to work with transverse electro- magnetic fields (the non-covariant 
electro- magnetic Dirac observables, the only ones which can be explicitly determined). 

The main aim of this paper is to write explicitly the Hamilton equations of ADM tetrad 
gravity coupled to this type of matter in the York canonical basis of Ref. [6] in the class of 
Schwinger time gauges, in which the tetrads are adapted to the 3+1 splitting of space-time 
(the time-like tetrad coincides with the unit normal to the 3-space S r and the three spatial 
tetrads have a conventional orientation). In these gauges the gravitational field is described 
by 10 configuration variables and 10 conjugate momenta like in canonical metric gravity. 
Since there are still 8 first-class constraints there are 8 gauge variables (inertial effects): 7 
configuration variables (the lapse and shift functions and 3 variables describing the freedom 
in the choice of the 3-coordinates in E T ) and a momentum variable (the trace 3 K(t, a r ) of 
the extrinsic curvature of the 3-space E r , describing the freedom in the clock synchronization 
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convention). One configuration variable (the element of 3- volume in the instantaneous 3- 
space E T , canonically conjugate to 3 K) and 7 momenta are determined by the 8 constraints. 
The remaining two pairs of conjugate variables describe the physical tidal effects contained 
in the gravitational field (they are only 3-scalars and partial Dirac observables with respect 
to 10 of the 14 constraints of tetrad gravity). 

As we shall see, the extrinsic curvature 3 K rs {r, a r ) of the instantaneous 3-spaces S T , 
viewed as sub-manifolds of the space-time, depends upon the matter, the tidal variables, the 
three inertial gauge variables for the choice of the 3-coordinates and, finally, upon the inertial 
gauge variable 3 K{r, a r ) (its trace). This last variable, the so-called York time, is what 
remains of the special-relativistic gauge freedom in the choice of the clock synchronization 
convention in general relativity. Once these gauge freedoms are fixed the extrinsic curvature 
and, as a consequence, the instantaneous 3-spaces are dynamically determined for every 
solution of Einstein equations with the allowed boundary conditions and with Cauchy data 
compatible with the constraints. 

In canonical ADM tetrad gravity the Dirac Hamiltonian is equal to the weak ADM energy 
Eadm (its form as a volume integral over the 3-space S T ) plus constraints. As a consequence, 
there is not a frozen picture like in spatially compact without boundary space-times. We 
will see that the gauge momentum (the inertial York time 3 K(r,a r ); its time nature is a 
reflex of Lorentz signature) gives rise to a negative kinetic term in Eadm vanishing only in 
the gauges 3 K(r, a r ) = 0. 

In Schwinger time gauges the ADM Hamilton equations in the York canonical basis are 
naturally separated in four groups. Besides the Hamilton equations for matter and for the 
tidal variables, there are four contracted Bianchi identities (implying the r-preservation of 
the super-Hamiltonian and super-momentum constraints) and the Hamilton equations for 
the four inertial gauge variables describing the 3-coordinates in E r and the clock synchroniza- 
tion convention (if we add four gauge fixings for these variables, these Hamilton equations 
determine the shift and lapse functions). 

Then we will define a family of non-harmonic gauges, the 3-orthogonal ones, in which the 
3- metric in the 3-spaces S T is diagonal but the inertial gauge variable 3 K(r,a r ), fixing the 
clock synchronization convention, is equal to an arbitrary numerical function F{r,a r ). The 
restriction of the Hamilton equations to the 3-orthogonal gauges is given explicitly. 

It is in this family of gauges that we will define a linearization of ADM canonical tetrad 
gravity plus matter in the second paper [19], to obtain a formulation of Hamiltonian post- 
Minkowskian gravity (without post-Newtonian expansions) with non-flat Riemannian 3- 
spaces and asymptotic Minkowski background. In particular this will allow to study the 
dynamical effects of the inertial potential 3 K(t, a r ) and to see whether it can, at least 
partially, mimic the effects attributed to dark matter. 

In Section II there is a review of ADM tetrad gravity, of the York canonical basis with 
new results beyond Ref.[6], showing the relevance of the Eulerian observers associated to a 
3+1 splitting of space-time, and of the ADM Poincare' charges. 

In Section III we give the action of our system in Subsection A and we evaluate the 
Dirac Hamiltonian and the constraints of the Hamiltonian formulation. In Subsection B 
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we give the Hamilton equations of the particles. In Subsection C, after having evaluated 
the Hamilton equations of the electro-magnetic field, we define the non-covariant electro- 
magnetic radiation gauge and give the restriction to it of the Hamilton equations of the 
particles and of the transverse electro-magnetic field. In Subsections D and E we give the 
Dirac Hamiltonian, the constraints and the ADM Poincare' charges in the York canonical 
basis, while in Subsection F we define the Schwinger time gauges. Then in Section IV we 
give the Hamilton equations in the York canonical basis in these gauges. 

In Section V we discuss the formulation in the York canonical basis and in the Schwinger 
time gauges of some of the most used gauges for canonical gravity, included the Hamiltonian 
formulation of the harmonic gauges for Einstein equations. 

In Section VI we give the restriction of the Dirac Hamiltonian, of the constraints and of 
the Hamilton equations to the above defined family of 3-orthogonal Schwinger time gauges. 

In the Conclusions there are some final remarks. 

In Appendix A there is a discussion of the contracted Bianchi identities and a comparison 
with the standard ADM Hamilton equations. 

In Appendix B there are the calculations needed for the Hamilton equations in the 
Schwinger time gauges, while in Appendix C these calculations are restricted to the 3- 
orthogonal gauges. 
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II. REVIEW OF TETRAD GRAVITY AND OF THE YORK MAP 



A. Tetrads and Cotetrads 

We use radar 4-coordinates a A = (a T = T;a r ), A = r,r, adapted to an admissible 3+1 
splitting of the space-time and centered on an arbitrary time-like observer: they define a non- 
inertial frame centered on the observer, so that they are observer and frame- dependent. The 
instantaneous 3-spaces identified by this convention for clock synchronization are denoted 
E T . 

The 4-metric a Qab has signature e (H ) with e = ± (the particle physics, e = +, 

and general relativity, e = — , conventions). Flat indices (a), a = o,a, are raised and 

lowered by the flat Minkowski metric 4 ^( Q )(/3) = e (H ). We define 4 ^( a )(6) = — e $(a)(b) 

with a positive-definite Euclidean 3-metric. From now on we shall denote the curvilinear 
3-coordinates a r with the notation a for the sake of simplicity. Usually the convention of 
sum over repeated indices is used, except when there are too many summations. 

We shall work with the tetrads a E^(t,ct) and the cotetrads 4 E A a \r, a). To rebuild 
the original tetrads 4 E^(r, a) = z A (T,d t ) 4 E A Y JT,d t ) we must know explicitly the embed- 
ding z^{r,a) of the instantaneous 3-spaces, so to be able to evaluate the transformation 
coefficients z a (t, a). 

General tetrads 4 E A t ~ ) (r, a) and cotetrads 4 E A a ^ (r, a) are connected to the tetrads and 

cotetrads adapted to the 3+1 splitting (the time-like tetrad is the unit normal l A to S r ) by 
a point-dependent standard Lorentz boost for time-like orbits acting on the flat indices 1 



o o \P) 

4 K) = 4 E W L^\ a) ( na) ), 4 E f = LW^W 'E a , 

, , («) o 03) 

4 g AB = 4 4 a) \ am 4 Ef = 4 E A \ am 4 E B , (2.1) 

where the last line gives the resolution of the 4-metric in terms of cotetrads. 

The adapted tetrads and cotetrads (in these so called Schwinger time gauges) have the 
expression (N(t, a) — 1 + n(r, a) > 0, with n(r, a) vanishing at spatial infinity (absence 



In each tangent plane to a point of S T the point-dependent standard Wigncr boost for time-like Poincare' 
orbits L(«\ 0) (V( Z (a)); £) = *;g+2eFW W ) ^ 

o(«) 

sends the unit future-pointing time-like vector V = (1;0) into the unit time-like vector V (a > = 
i E A a) l A = ^^/l + J2 a l -P\ a y^ a " 1 = — e( P(a)j, where l A is the unit future-pointing normal to E T . We 
have L _1 (<p( a )) = 4 ry L T (ip^) 4 i] = L(— <f( a ))- As a consequence, the fiat indices (a) of the adapted 
tetrads and cotetrads and of the triads and cotriads on E T transform as Wigner spin-1 indices under 
point-dependent SO(3) Wigner rotations Rr a \n,-\(V (z(a)); A(z(a)) ) associated with Lorentz transforma- 
tions A^^(z) in the tangent plane to the space-time in the given point of E T . Instead the index (o) of 
the adapted tetrads and cotetrads is a local Lorentz scalar index. 
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of supertranslations), so that N(t, a) dr is positive from E T to £ T +<i T , is the lapse function; 
N r (r, a) = n r {r, a), vanishing at spatial infinity (absence of super-translations), are the shift 
functions ) 



A E°a = (1 +n) (1;0) = el A , 4 E A = (n (a) ; 3 e (a)r ), (2.2) 

where 3 ej^ and 3 e( a ) r are triads and cotriads on S r and ri( a ) = n r 3 e r ^ = rf 3 e( a ) r 2 are 
adapted shift functions. 

The adapted tetrads 4 E^ are defined modulo SO (3) rotations 4 E^ = R^ b )(a^ e )) 4 E^, 
3 e r u) = R(a)(b)(oc( e )) 3 ^\b)' w here ot{ a ){j, a) are three point-dependent Euler angles. After 
having chosen an arbitrary point- dependent origin a( a )(r, a) = 0, we arrive at the following 
adapted tetrads and cotetrads [n^ a ) — Ylb n 0) R(b)(a)( a (e)) } 

4 E A } = 4 E° ] = (1 + n ) (1;0) = el A , 4^ = (n (a) ; 3 e (a)r ), (2.3) 

which we shall use as a reference standard. 
Then Eqs.(2.1), namely 

*Ef a) = % o) H°\ a) (<p {e) ) + % b) Rf b)(a) (a {c) ) LW (a) ( V(c) ), (2.4) 

show that every point-dependent Lorentz transformation A in the tangent planes may be 
parametrized with the (Wigner) boost parameters </?(„) and the Euler angles ct( a ), being the 
product A = RL of a rotation and a boost. 

The future-oriented unit normal to S r and the projector on S T are 



Since we use the positive-definite 3-metric #(«)(&)) we shall use only lower flat spatial indices. Therefore 
for the cotriads we use the notation 3 ef^ d = 3 e^ a y r with £( )(6) = 3e [ a ) 3e (b)r- 
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l A = e(l + n)(l;0), 4 g AB l A l B = e, 

lA = ^ = ik ( 1; ~ n 1 = TT-n ( 1; 

z h B = S *-el A l B , 3 K = 3 h; = 0, % = n {a) 3 e[ a)) 3 h r s = 5 r s , 

3* 37 3 37 33 

i^tt ^^(a)^(a); forr ^^(o) ^(a)r? "rs ^ ^(a)r ^(a)s? 

3 ^ T = 3 r = 0, 3 h rs = -e 3 e r {a) % y 

(2.5) 



B. The 4- metric and the Canonical Variables. 

The 4-metric has the following expression 

Vr = e[(l +n) 2 - 3 g rs n r n s ] = e [(1 + n) 2 - n (a) ], 
4 9rr = -en r = -fn( (l ) 3 e( fl ) n 

4 _ _ 3 
Qrs C (y'rs; 

flVs = Mrs — e(a)r e (a)s, 9 = ^ = ^U) e (a)> 



3 r 



4 tt c 4 rr 



e (a) ™(a) 



(1 + n) 2 ' y (1 + n) 2 (1 + n) 2 

^/3„rs n 11 . ^3 r 3„s /r "(«/""') , 3 - 

^ = " e (2 ~ (1 +w ) 2 ) = ~ e e W e W (<>(«)(&) ~ (1 +n y h 9 =h, 



4^rs r (3„rs 71 71 \_ , 3r 3„s ^ n ( a ) n ( b )\ 3 „rs _ h rs 

. . .. . g (a) . ... ... 

V r V s - W^-tt^, 3 ^ = 7 = ( 3 e) 2 , 3 e = cfet 3 e (a)r , 

(1 + n) z 



= = = V7 (1 + n) = 3 e (1 + n). (2.6) 

The 3-metric 3 g rs = h rs has signature (+ + +), so that we may put all the flat 3-indices 
down. We have 3 g ru3 g us = 5 r s (h ru h us = 5 r s ), d A 3 g rs = - 3 g ru3 g sv d A 3 g uv . 

The conditions for having an admissible 3+1 splitting of space-time are: 

a) 1 + n(r, a) > everywhere (the instantaneous 3-spaces never intersect each other); 

b) the M0ller conditions [10, 16], which imply 
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i) e 4 g TT > 0, i.e. (1 + n) 2 > 3 g rs n r n s (the rotational velocity never exceeds the velocity 
of light c, so that the coordinate singularity of the rotating disk named " horizon problem" 
is avoided); 

ii) e 4 g rr = — 3 g rr < (satisfied by the signature of 3 g rs = h rs ), 4 grr 4 9ss — ( 4 fiVs) 2 > 
and dete 4 g rs = —det 3 g rs < (satisfied by the signature of 3 g rs ) so that det A gAB < 0; 
these conditions imply that 3 g rs = h rs has three definite positive eigenvalues A r = A 2 in the 
non-degenerate case without Killing symmetries, the only one we consider; 

c) the space-time is asymptotically Minkowskian with the instantaneous 3-spaces orthog- 
onal to the ADM 4-momentum at spatial infinity: they are non-inertial rest frames of the 
3-universe (isolated system), there is an asymptotic Minkowski background 4- metric and 
there are asymptotic inertial observers whose spatial axes e% are identified by the fixed 
stars. 

As said in Ref. [6], in ADM canonical tetrad gravity the 16 configuration variables are: the 
3 boost variables ip( a )( T , &)\ the lapse and shift functions n(r, a) and 71( )(t, a); the cotriads 
3 e (a)r (r, a). Their conjugate momenta are 7^ (t, ct) , 7r n (r, a), 7r n(o) (r, a), 3 7r[ a) (r, ct). There 
are 14 first-class constraints: A) the 10 primary constraints ^ (r, a) 0, n n (r, a) ps 0, 
7r n(a) (r, a) ~ and the 3 rotation constraints M( a ) (r, a) ~ implying the gauge nature of 
the 3 Euler angles a( )(r, a); B) the 4 secondary super-Hamiltonian and super-momentum 
constraints H(r, a) 0, 7i( a )(r, a) « 0. As a consequence there are 14 gauge variables (the 
inertial effects) and two pairs of canonically conjugate physical degrees of freedom (the tidal 
effects). 

At this stage the basis of canonical variables for this formulation of tetrad gravity, natu- 
rally adapted to 7 of the 14 first-class constraints, is 



<P(a) 


n 


n {a ) 


3 e( a )r 1 




w 


7T„ W 


^"(a) 


=3 


K> 1 



(2.7) 



From Eqs.(5.5) of Ref. [3] we assume the following (direction-independent, so to kill 
super-translations) boundary conditions at spatial infinity (r = y/^2~ r (o" r ) 2 ; e > 0; 
M = const.): n(r,a) ^ r ->oo 0(r~ (2+e) ), 7r n (r, a) ^r^oo 0(r~ 3 ), n (a) (r, a) ->>^oo 
0(r" e ), vr„ (a) (r,CT) ^ r ^oo 0(r- 3 ), </? (a) (T,<7) ^ r ^oo 0(r-( 1+e )), 7r V(o) (r,5) ^ r ^oo 0(r- 2 ), 

3 e (a)r (T,,7)^oc (l + %)8ar + 0(r-W),**l a) (T,#) 0(r- 5 / 2 ). 
C. The York Canonical Basis. 

In Ref. [6] we studied the following point canonical transformation (it is a Shanmugad- 
hasan canonical transformation [20]) on the canonical variables (2.7), implementing the York 
map of Refs.[17, 18] and identifying a canonical basis adapted to the 10 primary first-class 
constraints . It is realized in two steps and leads to the following York canonical basis 
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3 e(a)r 1 
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¥>(a) "(a) 
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vr n « 7r S(a) « 


7T?> 




n„| 



(2.8) 

where 3 e( a ) r = ^ 6 3 e( b ) r R( b )(a)(oc( e )) (with conjugate momenta 3 7T( a) ), n^a) = 
^ b ri(b) R(b)( a )((X(e)) are the cotriads and the shift functions at a( a )(r, a) = after the ex- 
traction of the rotation matrix R^ a ^(a^(r, <?)), see after Eq.(2.2). 

With the first canonical transformation we extract the 3 angles a^(r, a) from the cotriads 
and we Abelianize the rotation constraints 3 M( a )(r, a) ~ 0, replacing them with the momenta 

conjugate to the angles, n^j = —J2 b 3 M^ A( b )(a) (at( e )) ~ 0. The 0(3) Lie algebra-valued 
Cartan matrix Ai m (9 n ) [B = A" 1 } is defined in Ref.[3]. 

The second canonical transformation is based on the fact that the 3-metric 3 g rs is a 
real symmetric 3x3 matrix, which may be diagonalized with an orthogonal matrix V(9 r ), 
V^ 1 = V T (^2 U V ua V u b = S a b, V ua V va = 5 UV , ^2 UV e wuv V ua V v b = e a bc V cw ), det V = 1, 
depending on 3 Euler angles 9 r 3 4 . The gauge Euler angles 9 r give a description of the 3- 
coordinate systems on S T from a local point of view, because they give the orientation of 
the tangents to the three 3-coordinate lines through each point (their conjugate momenta 
are determined by the super- momentum constraints). 



In the York canonical basis we have (from now on we will use V ra for V ra (9 n ) to simplify 
the notation) 



3 Due to the positive signature of the 3-metric, we define the matrix V with the following indices: V ru . 
Since the choice of Shanmugadhasan canonical bases breaks manifest covariance, we will use the notation 
Vua = J2 V Vuv 8 v (a) instead of V u ( a y We use the following types of indices: a = 1, 2, 3 and 5 = 1,2. 

4 A similar diagonalization of the 3-mctric has been considered also in Rcf.[21] for numerical grav- 
ity purposes. Instead of the Euler angles 9 n one uses an exponential representation of the ro- 
tation matrices V ra (0 n ) = e^™ tTam (— oo < Ci j C2 3 C3 < so that one gets 3 g rs = 
J2 a eX™ *r.mCm(fl') e 0a e E„ UanMO') with e K = ^ Q\. By defining 77 fe = ( k /( with ( = ({ + C| + C3 
(—1 < r/k < 1), one gets V ra — S a b cos C + X) c e °-bc Vc sin ( + i] a r)b (1 — cos £) and a similar expression for 
3 g rs . Then d 3 g rs is expressed in terms of d(f> a and d( a and, when the 3 eigenvalues are distinct, one can 
get d<p a and d( a in terms of d 3 g rs . 
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9rr = € [(1 + n) 2 - U U 

a 

9rr = -e E n (a) 3 e (a)r = -e0 1/3 ^ V ™, 

a a 

l 9rs = -e 3 9rs = -tJ2 V ruKS uv V v T s = -€ ^ (V m A°) (V sa A a ) = 

«u a 

= ~ e E 3 e(a)r 3 e(a) s = -e 3 e (a)r 3 e (a)s = -e 4 3 # rs = -e 2/3 Q 2 V ro V« 

a a a 

Aa = 6 «> = 2 = 1/3 Qa, Qa = e^' ^ = (?*\ 



<t> = 6 = Vl = Vdet 3 g = 3 e = ^Ai A 2 A 3 = Ai A 2 A 



3- 



3 e( a )r = ^ -R(a)(6)("(e)) 3 e(fe) r , 3 e (a ) r = 1/3 Q a Ka- 

M«) = E ^ ^(a) = ^^Q- 1 K, 

Ma) = E ^©("W)^)' 

^\a) « r 1/3 [V™ (0 ^ + E TSa + 



_l n- 1 Yll €alt Vwt r -(*) 



^ = - > y A m i[0 )e mir e (a ); 7r (a) , 



imra 

" ~~ \2k(t 3 ■''( 

rab 



Ha = E 7a6 3 ^ a) ^(a)(6)(tt(e)) 3 e {b ) r . (2.9) 



rab 



The set of numerical parameters 7 So satisfies [1] £) u 7 Su = 0, Y. u 7a« 7S« = *a5, E a 7a« 7a« = 
5 U „ — |. Each solution of these equations defines a different York canonical basis. 

We only consider 3-metrics with 3 distinct positive eigenvalues \ r = A 2 : the degenerate 
cases should be treated by adding by hand the constraints A a — A5 pa or Ax pa A 2 ~ A 3 
and by studying the resulting constraint algebra. 

The assumed boundary conditions given after Eqs.(2.2) imply A (r, a) = 
(4> 1/3 Qa)(r,a) ^™ 1 + f + ^_ + o(r~ 3 ) and 0(r, a) 1 + 0{r~ l ). Moreover 
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we must have nf\r, a) ^ r ^oo 0(r 4 ), since the requirement A a (r, a) ^ A(,(r, <?) for a ^ b, 
needed to avoid singularities, implies a a ^ a b for a ^ b in their asymptotic behavior, so that 

we get ( r '^) = (QbQ^-QaQb 1 ) (t,#) ^r^oc 2 {a b -a a y As a consequence, 

we have 7T/°x (r, (?) — > r - >-oo 0(r~ 5 / 2 ). Also the angles a( a )(r, a) and 0'(t, (?) must tend to zero 
in a direction-independent way at spatial infinity. We also have 7r^(r, a) ^ r ^oo 0(r~ 5 / 2 ) 
at spatial infinity. 

In Eq.(2.9) the quantity 3 K(r, <?) is the trace of the extrinsic curvature 3 K t . s (t, (?) of the 
instantaneous 3-spaces E T . In the York canonical basis the extrinsic curvature 3 K rs , the 3- 
spin connection 3 u> r ( a ) = 3 coV(a)U (e) =o ( see Eq.(B17) of Ref. [6]) and the 3-Christoffel symbols 
have the following expression [6] 



c 3 3 e 



ATT G \ I" ^3 _ 3 _ 3_ 3 _ \ 3 _ _ u 3 _ 3 _ 3 _ _ u 

2 , I l e (a)r e(b) s + e( ) s e(6) r I e( a ) u 7T( b )- e( a ) r e( ) s e(6)„7T( 6 ) 

r 1/3 ( E ^ ^» ^~ [ 2 E r 6a n s - 4> + 



afeu 

4ttG ~ 



a 6 



QbQ^-QaQt 



■ Co' 



r(o) 



Ha)(6)(c) ^r(6)(c) 



fee 



<9r 3 e( c ) u - <9 U 3 e( c ) r + E 3e {"c) 3e (d)r d v 3 e( d ) u 

dv 

= E e ( a )( fe )( c ) ^ _ ^ Q° ^ ^ ^ + E %c ^) + 

bcu a 

+ i Q, 1 Q c (d r V uc - d u K c ) + i E ^ ^ ^ ^ ^ ^ 



abs 

+ E 75a ^™ (^ Ma 9v R& + Vva 9u Rd ) ~ E Tbb Qa 2 V ra V sa V ub V vb 8 S R b + 

aa l abs 

+ \ E Vra (^ v ™ + 9 « Ka ) + 



+ « E 2 ^ 2 ^« ^< 



afes 



Kfe dv V sb - d s V vb ) +V vb [ d u V sb - d s V x 



'ub 



E 3 ri = r'dui 



(2.10) 
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The previous sequence of canonical transformations realizes a York map because the 
gauge variable ir^ (describing the freedom in the choice of the instantaneous 3-spaces S T ) 
is proportional to York internal extrinsic time 3 K. It is the only gauge variable among the 
momenta: this is a reflex of the Lorentz signature of space-time, because and 6 n can be 
used as a set of 4-coordinates [7]. 

Its conjugate variable, to be determined by the super-hamiltonian constraint, is <fi = 6 = 
3 e, which is proportional to Misner's internal intrinsic time; moreover is the 3-volume 
density on S T : Vr = J R d 3 a<f) 6 , R C S T . Since we have 3 g rs = 4> 2 ^ 33 g rs with det 3 g rs — 1, 
is also called the conformal factor of the 3-metric. 

The two pairs of canonical variables R a , IT^, a = 1,2, describe the generalized tidal 
effects, namely the independent degrees of freedom of the gravitational field. In particular 
the configuration tidal variables R a depend only on the eigenvalues of the 3-metric 5 . They 
are Dirac observables only with respect to the gauge transformations generated by 10 of the 
14 first class constraints. Let us remark that, if we fix completely the gauge and we go to 
Dirac brackets, then the only surviving dynamical variables Ra and U a become two pairs of 
non canonical Dirac observables for that gauge: the two pairs of canonical Dirac observables 
have to be found as a Darboux basis of the copy of the reduced phase space identified by 
the gauge and they will be (in general non-local) functionals of the R a , U a variables. 

~ (0) 

Since the variables <p and ir\ are determined by the super-Hamiltonian and super- 
momentum constraints, the arbitrary gauge variables are a^ a ), (p( a ), @\ ^ n an d n( a )- As 
shown in Refs.[5], they describe the following generalized inertial effects: 

a) a( a )(r, a) and (p( a )(j, a) are the 6 configuration variables parametrizing the 0(3,1) 
gauge freedom in the choice of the tetrads in the tangent plane to each point of S T and 
describe the arbitrariness in the choice of a tetrad to be associated to a time-like observer, 
whose world-line goes through the point (r, a). They fix the unit ^-velocity of the observer 
and the conventions for the orientation of gyroscopes and their transport along the world-line 
of the observer. 

b) 9 i (r, a) [depending only on the 3-metric, as shown in Eq.(2.9)] describe the arbitrariness 
in the choice of the 3-coordinates in the instantaneous 3-spaces S T of the chosen non-inertial 
frame centered on an arbitrary time-like observer. Their choice will induce a pattern of 
relativistic inertial forces for the gravitational field , whose potentials are the functions 
V ra (0 l ) present in the weak ADM energy Eadm given in Eqs.(3.14). 



5 If we consider the eigenvalue equation for the 3-metric 3 g rs of determinant one, we identify the following 
two 3-scalars depending only on the tidal variables R a ' i) Tr 3 g rs = ^2 r 3 g rr = Yl a Qa (the- sum of the 
eigenvalues); ii) 3 g n 3 g 22 - 3 g\ 2 + 3 g 2 2 3 533 - 3 & + 3 .933 3 .9n - 3 gli =QlQt + Ql Q\ + Q\ Q\ (the sum 
of the possible products of two eigenvalues). 

This suggest the possibility of a point canonical transformation from i? a , II a to 3-scalar tidal variables 
X{R- a ), Il x , Y{R- a ), U Y with 



711 #I+72i #5 



X = e L'" ' ' 'J +e 

(711+712) -Rl + (75l+722) R2 



712 #1+722 #5 



71 3 #1+723 #2 



Tr 3 g rs = 4>- 2/3 Tr ' 



'fh 



Y = e 
= ^- 4 /3 



+ e 



(712+713) #l + (722+723) #5 



(7i3+7ii) #I+(723+75l) #5 



3 322 " 3 Z?12 + 3 522 3 333 " 3 <?! 3 + 3 533 3 ffll ~ 3 <?31 
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c) fj( )(r, a), the shift functions appearing in the Dirac Hamiltonian, describe which points 
on different instantaneous 3-spaces have the same numerical value of the 3-coordinates. They 
are the inertial potentials describing the effects of the non-vanishing off-diagonal components 
A g Tr {r,a) of the 4- metric, namely they are the gravito -magnetic potentials 6 responsible of 
effects like the dragging of inertial frames (Lens-Thirring effect) [22] in the post-Newtonian 
approximation. The shift functions are determined by the r-preservation of the gauge fixings 
determining the gauge variables 6 % {r, a). 

d) 7r^(r, a), i.e. the York time 3 K(r,a), describes the non-dynamical arbitrariness in 
the choice of the convention for the synchronization of distant clocks which remains in the 
transition from special to general relativity. As said in the Introduction, the choice of the 
shape of the instantaneous 3-space as a sub-manifold of space-time (a pure gauge choice in 
special relativity) is dynamically determined by the chosen solution of Einstein's equations 
after the fixation of the gauge variables 3 K(r, a) and 9 l (r, a). Since the York time is present 
in the Dirac Hamiltonian 7 , it is a new inertial potential connected to the problem of the 
relativistic freedom in the choice of the instantaneous 3-space, which has no non-relativistic 
analogue (in Galilei space-time time is absolute and there is an absolute notion of Euclidean 
3-space). Its effects are completely unexplored. 

e) n(r, a), the lapse function appearing in the Dirac Hamiltonian, describes the arbitrari- 
ness in the choice of the unit of proper time in each point of the simultaneity surfaces £ r , 
namely how these surfaces are packed in the 3+1 splitting. The lapse function is determined 
by the r-preservation of the gauge fixing for the gauge variable 3 K(t, a). 

The gauge variables 6 1 (t,(t), n(r,cx), n^(r,a) describe inertial effects, which are the 
the relativistic counterpart of the non-relativistic ones (the centrifugal, Coriolis,... forces 
in Newton mechanics in accelerated frames) and which are present also in the non-inertial 
frames of Minkowski space-time [10]. 

D. The Expansion and the Shear of the Eulerian Observers. 

Let us now consider the geometrical interpretation of the extrinsic curvature 3 K rs of the 
instantaneous 3-spaces £ r in terms of the properties of the surface-forming (i.e. irrotational) 
congruence of Eulerian (non geodesic) time-like observers, whose world-lines have the tangent 
unit 4- velocity equal to the unit normal orthogonal to the instantaneous 3-spaces S r . If we 
use radar 4-coordinates, the covariant unit normal el a — (1 + n) (1; 0) of Eqs.(2.2) has the 
following covariant derivative (see for instance Ref. [23]) 



6 In the post-Newtonian approximation in harmonic gauges they are the counterpart of the electro-magnetic 
vector potentials describing magnetic fields [22]: A) N = 1 + n, n d = — |# $g with <&g the gravito- electric 
potential; B) n r d = |# Aq t with Aq t the gravito-magnetic potential; C) Eq t = d r 3>g — d T (| Ao r ) (the 
gravito- electric field) and Bq t = e ruv d u A Gv — cVLq t (the gravito-magnetic field). Let us remark that in 
arbitrary gauges the analogy with electro-magnetism breaks down. 

7 See Eqs.(3.43) and (3.44) for its presence in the supcr-Hamiltonian constraint and in the weak ADM 
energy, and Eqs.(3.41) for its presence in the super-momentum constraints. 
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4 Va e Ib — £ I a 3fl B + &AB + ^ h AB — ^ab — e I a 3 °b + 3 K AB , 



} K AB = 3 K rs b\ b s B , b r A = 5 r A + 3 e( a ) 5 A , h AB = g AB — e l A l B . 



[2.11) 



The quantities appearing in Eqs.(2.11) are: 

a) the acceleration of the Eulerian observers 

3 n A — / B4 Y7 lA _ 4 n AB3 n 3 n _3 n ir 

a —l V B l> — 9 OjB, CLA — drOji, 

3 a r = d r In (1 + n), 3 a T = 3 a r 3 e[ a) 3 e u {a) n u = 3 a r 3 e\ a) n (a) , 

3 a r = -e 3 e r (a) 3 el a) 3 a s , V = 0; (2.12) 

b) the vorticity or twist (a measure of the rotation of the nearby world-lines infmitesimally 
surrounding the given one), which is vanishing because the congruence is surface-forming 

u AB = -uj ba = e - (l A 3 a B - l B 3 a A ) - | ( 4 Va Ib - 4 V B I a) = 0, 



co AB l B = 0, cu A = ^- 4 n ABCD cu BC l D = 0; (2.13) 

c) the expansion 8 , which coincides with the York external time, is proportional to the 
Hubble parameter H 9 and to the dimensionless (cosmological) deceleration parameter q = 
3l A4 V A l-l = -39- 2 l A d A 9-l, 

f) 4 X7 1 A , 3 K 4n ° H[a)r ^ , 127T °ir 
6= V A l =-eK = = - e _^_7r-, 

H = ±6=jl AA V A l = -e^-^, g = 3/ A4 V A i-l; (2.14) 



It measures the average expansion of the infmitesimally nearby world-lines surrounding a given world- line 
in the congruence. 

I is a representative length along the integral curves of 4 E^ , describing the volume expansion (contraction) 
behavior of the congruence. 
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d) the shear 
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(TAB = °BA = - ^ (V IB + 3 dB I A) + ^ (*V A h + ^ B U) ~ g 9 3 h AB = 

= ( 3 K rs - 1 3 g rs 3 K) b A b B , A g AB a AB = 0, a AB l B = 0. (2.15) 

By explicit calculation we get the following components of the shear along the tetrads 
(2.2) 

o(«) oW) 

0AB = CT(q)(/3) -^B = fi'AC* fi'BD & , 

Ott = n(a) n(b) a (a)(b) = ( 3 K rs - - 3 g rs 3 K) n (a) 3 e r ^ 3 e s ^, 



' TV 
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a (a)(b) 3 6(a)r 3 e(b) s = 3 if rs ~ ^ 3 <7rs 3 -^, 



a = a = U, cr = e (a) e (b) cr (a)(b) , 
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<7(o)(a) 


= &AB 




4 ^(a) = 






o A 


o B 


(7(a)(6) 


= ^AB 


4 E (a) 


4 ^(b) = 



(1 + w )2 [^tt - 2 n (a) 3 e( a) a rr + n (a) n (6) 3 e[ a) 3 e( b) <r„], 
(cr rs - 3 e( fe) cr rs ) 3 e( a) , 



3 3 -s 



4 

°"(o)(o) = 0, 
°"(o)(a) = 0, 

^( a)(6 ) = a ma) = { 3 K rs - i 3 g rs 3 K) 3 e\ a) 3 e s (b) , £ <7 (a)(o) = 0. (2.16) 
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It measures how an initial sphere in the tangent space to the given world-line, which is Lie-transported 
along the world-line tangent (i.e. it has zero Lie derivative with respect to d^), is distorted towards 
an ellipsoid with principal axes given by the eigenvectors of cr%, with rate given by the eigenvalues of 
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cT( a )(b) depends upon 6 r , 0, R a , and n s . 

As a consequence we have 11 (n r \ s is the covariant derivative on S T ; = means evaluated 
by using the equations of motion) 



K rs = — - 3 g rs 9 + (7(a)(6) 3 e( a ) r 3 e({,) s = 2 ^ 3 E ^ — - 05 ab + C(a)(&)) Qa <3& Ka Vsb, 

aft 

«9 T 3 9rs (2 = 9) 2 / 3 E q 2 [2 (1 0- 1 9 r + 9 r r«) v; a + d T (y ra k«) 

a 

= n r \ s + n s \ r - 2 (1 + n) 2/3 E Qa Q& Ka Kb ( - I 6> 5 afe + (7(a)(6)) ■ 



(2.17) 



By using Eqs.(2.9) and (2.10) we get 
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87TG ^ C 3 ~ sr^ 

'(7(a)(a) = — 2_^i laa^i-a, ~> Us = ~~^~Q ( P 7aaC(a)(a), 



„ I de/ 7 , 87TG CabtVwt ST^ D (0) 

\a)(b)Wb = <PV(a)(b)\a^b = 2^ H H~l - H H"! L B ™ K i ' 

<- tw V& V a Va Vfc i 

Ctab S( a )(6)|a^b, 
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(0) c 



8irG 

wtab 



4vr G a ^ 



(7(a)(6) QaQbV ra Vsb - 

a ab 
J2TaaTlaQlV ra V sa \, 



8ttG ~ x 



C 3 



2 de/ 1 1 /87rG\ 2 ~ 2 2 

CT = 2 ^ a(a)(b) f7(a)(b) = 2 l^ 3 " J I ^ s + 

afe a 



KulKu'l Kt)3K)'3 



, 2 ( wl w>i i Vw2 Vw ' 2 i_ 

H VQ2Q, 1 -Q,Q- 2 1 ) 2 {Q3Q1 1 - Q1Q3 1 ) 2 (QiQ^-QzQl 1 ) 2 - 

WW 

E^^-'^^r]- ( 2 - 18 ) 



11 The 3-scalars associated to the symmetric matrix 3 K rs are I = 3 K = —e6, II = det 3 K rs , III = 
3 Kn 3 K 22 - 3 K\ 2 + 3 K 22 3 K 33 - 3 Kl 3 + 3 K 33 3 K X1 - 3 K\ X . If 3 K rs = 3 K rs - | 3 g rs 3 K is the traceless 
extrinsic curvature, the 3-scalars // and 77/ may be replaced by II = del 3 K rs and III = 3 Kn 3 K 22 — 



3 K 2 12 + 3 K 22 3 K 33 - 3 K| 3 + 3 K 33 3 K U - 3 Kl 1 . 

12 



The canonical tidal variables R a , n a are 3-scalars, since they can be replaced with a canonical basis built 
with the 3-scalars X(R n ), Y(R a ) (see footnote 5) and II a = -g^j 4> J2 a 7aa <?(a)(a)- 
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Therefore the Eulerian observers associated to the 3+1 splitting of space-time allow a 
physical interpretation of some of the variables of the York canonical basis: 

1) their expansion 9 is the gauge variable determining the non-dynamical part of the 
shape of the instantaneous 3-spaces S r ; 

2) the diagonal elements of their shear describe the tidal momenta while the non- 
diagonal elements are connected to the variables nf\ determined by the super-momentum 
constraints. 



E. The Asymptotic ADM Poincare' Algebra, the Rest-Frame Conditions and the 
Center of Mass 

As explained in Ref.[l], following suggestions of Dirac in Ref. [24], the limit of the em- 
bedding z^ij, a) for the non-inertial rest-frame instant form of metric and tetrad gravity at 
spatial infinity in asymptotically Minkowskian space-times has the form 

z»(r,a) -^ co) (r) + 6f oo)r ^, (2.19) 

where (r) is an asymptotic inertial observer and = 6^ r and 6^ r are an asymptotic 
tetrad (denoted in the Introduction). 

The generators of the asymptotic ADM Poincare' group are [1] 



P(oo) ~ °(oo)A\ T ) ADM ~ '(oo) ADM + ± ADMi 

TV I* v _ v u , nfJ-v 

J (oo) — X (oo)^(oo) (oo) .r (oo) (oo) ' 



°(oo) — (oo) A (oo) B J ADM ~ 

~ itfoo) b (oo)r - ^(00) & (oo)r] J ADM + [ & (oo)r &(oo)s ~~ &(oo)r ^(oo)J '^A & DMi (2.20) 

where is the variable canonically conjugate to x^j. Padmi ^adm are the s ^ ron 9 (surface 
integrals at spatial infinity) ADM Poincare' charges in adapted radar coordinates: they are 
weakly equal to the weak (volume integrals over E T ) ADM Poincare' charges Padmi Jadm- 

Since is orthogonal to the asymptotic hyper-plane (2.19) due to the requirement 
of absence of super-translations [1], we can make the identifications = b^y. = 

P(oo)/ \J e P%o) = = (V 1 + h 2 ; h) and b^ r = e%(h), where and the space-like 4-vectors 

€%(h), orthogonal to are the column of the standard Wigner boost sending P M = Mch^ 
to its rest-frame form Mc (1; 0) (see Ref. [8] for this notation). 

Then for consistency the first of Eqs.(2.20) implies the rest-frame conditions 

Padm ~ Padm ~ 0) (2-21) 

and Mc « P\dm = E ADM /c. 
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From Eqs.(25) of Ref.[3] the weak or volume form of the ADM Poincare charges appearing 
in Eqs.(2.20) is 



Padm — ~ Eadm — I d 3 a 



+ 
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J ADM 



= J d 3 a [(a 



_s 3-pr \ 3ttud 



-(cr 5 -a g )M U (T,a). 



(2.22) 



They are the same weak Poincare charges of metric gravity, expressed in terms of cotriads 
3 e (a)r and their conjugate momenta 3 ir r {a) , by using 3 g rs = 3 e (a)r 3 e (a)s , 3 U rs = \[ 3 e r (a) 3 vr? } + 



3 e s (a) 3 ?rLJ (see Eq.(5.7) of Ref.fl]). The Christoffel symbols 3 r^ are built with the 3-metric 



'(a) "(a) 
9rs- 

In Eqs.(2.22) we use the notations A) Ai(T,a) for the internal matter mass density (in 
general metric-dependent); B) A / l r (r, a) for the internal matter momentum density (metric- 
independent and universal). See Section III, Eq.(3.9), for their explicit form: there we will 
give the modifications of the relevant formulas of Ref. [3] due to the presence of matter, 
because, in absence of derivative couplings to matter, they are independent from the type of 
matter: only the explicit form of the mass density M. (and also of the matter stress tensor 
T rs ) depends on the type of matter 13 . 

In Section III, Eq.(3.11), there is the form of the energy-momentum tensor T ab (t, a) for 
the type of matter considered in this paper. Since our formulation is equivalent to the ADM 



13 We have: 

A) For the Klein- Gordon field the Lagrangian is \ 4 g AB dAfds <p — m 2 tp 2 (r, (?) and we have M. (r, (?) = 

\ ?j- + 3 e ( 3e ( a ) 3e ( a ) dr^ds (f^j (t, (?) and M t (t, (?) = 7t(t, (?) d r ip(r, a), where 7r(r, <?) is the KG momen- 
tum. 

B) For perfect fluids [25] with Lagrangian coordinates a 1 (r, (?) and equation of state p = 
p(n,s) and fluid momenta ITi(r, a) we have A4(t, a) = 4> 6 X p((j)~ 6 X) + ([det (d r a 1 )] 2 



X 2 ) OP q* ] \ x =4>- 6 x <t) (X = X( 3 g rs , <f>, a 1 , IL) is the solution of a trasccndcntal equation depend- 
ing on the equation of state) and A4 r (r, a) = — Hi(r, (?) d r a l (r 7 a). For a dust we have A4(t, a) = 
\J M 2 c 2 [det (d r a* )} 2 + e ~ 2 Ea 7au fla n m d u a m Tl n d u a n {r, a). 
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one, and therefore to Einstein equations, we have A V aT ab (t,(j) = as a consequence of 
the Bianchi identities. 

In Ref. [26] it is said that the ADM Poincare' charges coincide with those arising from the 
Landau-Lifshitz energy-momentum pseudo-tensor for the gravitational field t^ L = t v ^ L [27] 
(see Ref. [28] for a detailed analysis). As a consequence the identities 4 V A T ab (t, a) = can 

be rewritten in the form of conservation laws d B 



(- A g) (t ab + tH ) 



o. 



In Ref. [10] we gave the form of the internal Poincare' generators of an isolated system in 
special relativity in both the inertial and non-inertial rest frames. By restricting Eqs.(2.20) 
to Minkowski space-time (a special solution of Einstein's equations when G = 0) and then to 
inertial rest-frames (where 4 gAB = 4 Vab and the embedding identifies the Wigner 3-spaces), 
it turns out that the ADM Poincare' generators become equal to the special relativistic 
ones for the given matter [8] in the inertial rest-frames. Also the non-inertial rest-frame 
Poincare' generators [10] could be recovered without the restriction to inertial rest-frames 
(now 4 g AB = A g AB [z]). 

As a consequence, this approach leads to the following visualization of the 3-universe 
(with its content of gravitational field and matter) contained in the instantaneous 3-spaces 
S r : it is an isolated system described by a decoupled external canonical non-covariant 
4-center of mass 5 M (r) (a non-observable decoupled pseudo-particle; it replaces the 

of Dirac proposal) carrying a pole-dipole structure with mass Mc ~ Eadm/c and a spin 
g r _ i e ruv J^" DM ~ 1 e ruv J^ DM . It is the same structure present in special relativity 
[8, 10] in the inertial and non-inertial rest-frame instant forms of dynamics. There, in 
non-inertial rest frames the embeddings tend at spatial infinity to the Wigner hyper-planes 
Z ^( T , a) = Y»(t) + e?(h) o r (the analogue of Eq.(2.19)), where Y*{t) = Y»{0) + r is the 
covariant non-canonical Fokker-Pryce center of inertia of the isolated system. Both x^(t) 
and Y^{t) are well defined functions (given in Ref.[8]) of r, Mc, S r , and of frozen (non- 
evolving) Jacobi data z, h, for the decoupled external center of mass (x^w = z/Mc is the 
non-covariant Newton- Wigner 3-position). 

5 M (r) carries also a universal external realization of the Poincare' algebra with generators 
PM = Mc/i", J ij = z i h j - z j K + e iju S u , K l = J oi = -y/l + h 2 z i + {gx . h) \ . In tetrad 

gravity these generators are Let us remark that the non-observability of this 

decoupled pseudo-particle is in accord with the viewpoint of Ref. [29] . 

Inside each instantaneous 3-space there is the isolated system and an internal realization 
of the Poincare' algebra, whose generators correspond to the weak ADM Poincare' charges. 



If, like in special relativity [8, 10], we eliminate the internal 3-center of mass inside the 
instantaneous 3-spaces S T with the following gauge fixings to the constraints (2.21) (together 
they form three pairs of second class constraints) 

J ADM ~ 0) (2.23) 

then we can choose the world-line of the non-canonical covariant Fokker-Pryce 4-center of 
inertia as origin of the 3-coordinates a r in the 3-spaces S T (instead of x 1 ?,). 
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Let us add a final remark. In Eqs. (12.2) and (12.3) of Ref.[l] and in Section 7 of Ref. [3] 
it is shown how to make a so-called Frauendiener- Sen- Witten transport (depending on the 
extrinsic curvature of E r ) of the asymptotic flat tetrads e^(h) to each point of E r : this 
allows to define a local dynamical compass of inertia A E^ to be compared in each point of 
E r with the tetrad A E^ a y These special tetrads have the following form (a( a ) are suitable 
Euler angles) 

4 K) = z "%) "> "cW, 3 e\ c) = R {c){d) (a {e) )% } - 5 a c , (2.24) 

with the special triads 3 e r ^ solution of the following Frauendiener equations [30] (a is a 
proportionality constant) 



3 V r % } = 3 V r % } = 0, 3 V r 3 e[ 3) = -a 3 K = -a ^, 

3 e(i) 3 ef 3) 3 V r 3 e (2)s + 3 e r {3) 3 ef 2) 3 V r 3 e (1)s + 3 e[ 2) 3 ef x) 3 V r 3 e {3)s = 0. (2.25) 

Therefore, these triads are formed by 3 vector fields with the properties: i) two vector 
fields are divergence free; ii) the third one has a non-vanishing divergence proportional to 
the trace of the extrinsic curvature of S T (on a maximal slicing hyper-surface, 3 K = 0, all 
three vectors would be divergence free); iii) the vectors satisfy a cyclic condition. In Ref. [30] 
it is shown these triads do not exist for compact S T and in general when there is nontrivial 
topology for E T . 

In conclusion, there are preferred ADM geometrical and dynamical Eulerian observers 
with unit 4- velocity 4 E^ and gyroscopes along the spatial axes 4 E^ a y The asymptotic 
world-lines of the congruence of these observers may replace the static concept of fixed stars 
in the study of the precessional effects of gravito-magnetism on gyroscopes (dragging of 
inertial frames) and seem to be naturally connected with the definition of post-Newtonian 
coordinates [31]. This congruence of time- like preferred observers is a non-Machian element 
of these noncompact space-times. 

These preferred tetrads correspond to the non-flat preferred observers of Bergmann [32]: 
they are a set of privileged observers (privileged tetrads adapted to the instantaneous 3- 
spaces E T ) of geometrical nature and not of static nature. Since they depend on the intrinsic 
and extrinsic geometry of E T , on the solutions of Einstein's equations they also acquire a 
dynamical nature depending on the configuration of the gravitational field itself. 



23 



III. THE ADM ACTION IN PRESENCE OF THE ELECTRO-MAGNETIC FIELD 
AND OF CHARGED SCALAR PARTICLES WITH GRASSMANN- VALUED 
ELECTRIC CHARGES AND SIGN OF THE ENERGY. 



Let us now describe N charged scalar particles and the electro-magnetic field coupled to 
the gravitational field in ADM tetrad gravity. 

A. The ADM Action and the Constraints 

The tetrad ADM action 14 for tetrad gravity (see Eq.(4) of Ref. [3]) plus the electro- 
magnetic field and N charged scalar particles with Grassmann-valued electric charges Qi and 
sign of the energy rji, depending on the configuration variables n{r, a), ri( )(r, a), ^( a )ij, a), 

3 e (a)r (r,a), A A (r,a), #(r), ^(t), 9f\r), is 15 
S Sg rav + S em + Sp ar t + SQ rassmann 

.3 



c 



16nG 



J drd 3 a{(l + n) 3 ee 



3 r 3s 3q i 
Wnl ^ L rs(c) ^ 



(a)(6)(c) C( a ) C(6) "rs(c) 



3 e 

+ 2(1 + n ) 1 )(«)C')(e)('0 3e wK)k ~~ , 9r 3 e(a)r) 3 e( d) (n (c) | s - <9 r 3 e (c) s )}(r, 5) - 



~ \ J drd 3 a 3 e(T,a) 



2 3 e^ 3 e? s 4 3 e!^ n fa ) 3 ef M 3 e^ M 

( a ) W p p _|_ (a) W (6) (6j j? j? , 



1 + n " 1 + n 

, 3 p r 3 p s 3 p « 3 p t, r/i ,_u r W( c ) 71(d) + g( c )(d) n (a ) n (6 ) 1 
+ e (a) e (fe) e (c) e (d) K 1 + "J <>( )(e,) d( c ) (d ) ^-j-^ \t ru t sv ^T j a)- 



14 Dimensions of the quantities appearing in this Section: [t] = [x* 1 ] = [a] = {rji] = [I], [Hi] = [P^] = [E/c] = 
M*" 1 ], [ 4 9 ] = ?g] = [n] = [„ (o) ] = [ 3 e (a)r ] = %] = [6,] = [0], [G = 6.7 1CT 8 cm 3 s~ 2 g^ 1 ] = [m" 1 lH~% 
[G/c 3 ] - [m" 1 *] « 2.510- 39 Sec /5, [5] = [ft] = [J AB ] = [mi 2 *" 1 ], [ 3 i?] = [ 3 fi rs(a) ] = [r 2 ], [V (a) ] = 
[ 3 X rs ] = [r 1 ], [ 3 ^ a) ] = [ 3 fl-] = [ml-H-^ [T AB ] = [M] = [M r ] = [H] = [H (a) ] = [m/" 2 *" 1 ]. 

15 The use of a positive 3-metric changes the sign definition of the particle momentum and the Poisson 
bracket sign with respect to Ref. [15]. G is the Newton constant. In Scrassmann we use the Planck 
constant h for dimensional reasons, since our regularization of the self-energies is considered a semi- 
classical approximation of a quantum theory, which exists for the electro-magnetic field (QED) but not 
yet for the gravitational field. 
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-^2 dTd 3 a5 3 {a,rji{T)) (m iC 
i=i J 

V* \/ ( 1 + n(r, a)) - ( 3 e (a)r (r, ct) ^ r (r) + n (a) (r, a)) ( 3 e (a)s (r, a) 77? (r) + n (a) (r, 



+ 



[A T (r,o : )+A r (r,a)77[(T)]) 



+ 



+ 



(3.1) 



While <5i is the Grassmann- valued electric charge of particle " i" , rji is its Grassmann- valued 
sign of the energy (particles with negative energy have the opposite electric charge: rjiQi). 



In the action S grav for the gravitational field 3 £l rs ( a 



(a) 



9 r 3 C^s(a) 



9 S 3 CoV( a ) 



e( a )(b)( c ) uo r (b) ^s(c) is the field strength associated with the 3-spin connection 



V(a) (T, CT) 



e ( „)( 6 )( c ) 



and ( 3 Go 1 )(a)(6)( c )(d) = #(a)( c ) + 5( a )(d)<fy>)(c) - 25( a )(&)5( c )(d) is the flat (with lower indices) 
inverse of the flat Wheeler-DeWitt super-metric 3 G ( a )^ c )( d ) = 5( a )(c) <fy>)(d) + £(a)(d) <fy>)(c) - 
5(a)(6) <J(c)(«o, 3 G , ( a )( 6 ) (e)(/ ) 3 G o( 1 e)(/)(c)(d) = 2 (<S (o)(c) 5 {b)(d) + 5 {a){d) 6 mc) ). 

J drd 3 a a/ - 4 g 4 g AC 4 g BD F AB F C d (t,v) has the previous 



The action S. 



_ _ 1 

em 4 

form because 4 g AC 4 g BD F AB F CD = 2( 4 g TT4 g rs - 4 g Tr 4 g TS ) F Tr F TS + A 4 g Tr 4 g su F TU F rs + 
y y r ru r sv 

The canonical momenta for the tetrad gravity variables and the electro-magnetic field are 
(we are in the canonical basis (2.7)) 



5S A 



DM 



5d T (p( a )(r, a) 



= 0, 



, ^ SSadm „ 
M r > <r) = — — — — = 0, 
do T n{r, a) 

8Sadm 
5d T n (a) (r, a) 



7T, 



0, 



3 < a )(r, s ) 



55 



ADM 



£<9 r 3 e (a ) r (T, a) 



167rG l 1 +n 



i G ^(^(bKcKd) 3 e(fe) 3e (d) ( n (c)\s ~ d T 3 e( c ) s )](r , 



cr 



8ttG 



[ e( K - e (c) e (c) A J e (o ) s J(r, crj, 
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o o T A t (t, a) 

r/ -\ & SaDM 



^ = (i^M.,M.,ft.-»'y)(^i. 



{A A (r, a),7r B (r, a')} = cr ] B 5 3 (a,a'), 

{n(r, <r), 7r n (r, a)} = 5 3 (a, a), 
{n {a) (T,a),7v n(b) (T,a)} = 5 {a){b) 5 3 (a,a), 
W(a){r,a),n nb) (r, a)} = S^S 3 ^, a), 
{ 3 e (a)r (r, a), 3 vr ( s b) (r, a)} = 5 {am 5 s r 5 3 (a, a), 

f3^r r „_ -?\ 3~s 



r( a) (r^')} = 3 e(r^) 3 e[ a) 
while the particle momenta are 



{ 3 e(r, a), 3 ^ a) (r, a')} = 3 e(r, a) 3 e[ a) (r, a) <5 3 (a, a'), (3.2) 



M T ) = 9L ^f^ d = ViKir(r) = A r (T,ffi(r)) + 
dr}l(T) c 

+ 7ft mi c 3 e {a)r (T,rf i (T)) ( 3 e (a)s (r, a) ?)|(r) + n (a) (r, a)) 

^(1 + n(r, a)) 2 - ( 3 e (a)r (r, a) t)[(t) + n (a) (r, a)) ( 3 e (a)s (r, a) t)?(t) +n (o) (r, a)) CT -* {r) ' 



) = J dipi dip* «i r (r), W(r), k^(t)} = 



4 



y m? c 2 + 3 e[ a) (r, a) Jk^j) - ^ A r (r, 3 ^ a) (r, a) (/c is (r) - ^ A s (r, a)) |^. (t) = 

mi [1 + n(r, a)] 

+ n(r, a)) 2 - ( 3 e (a)r (r, a) t)[(t) + n (a) (r, a)) ( 3 e {a)s {r, a) t)?(t) + n {a) (r, a)) 

(l + n) 3 ^ (n is {r)-9i A ^ 



\a=f)i(T) ? 



^jmjc^ + ^ 3efa (^(r) - % Z) (^(r) - f Z) 
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To define Ki r from k ir we have used the following property of Grassmann variables 16 

J d9i d9* = 0, J d9i d6* 9* 8i = 1. (3.4) 

The momenta of the Grassmann variables, implying second-class constraints, are 

*4 Q) W = -jgm = --2 e * (r)j \ (Q) * (r) = ^p = ~Y 9i (r)j 

{9f\r),n e Mr)} = R (Q) *(r), 7r flW ,.(r)} = 

3 J 

'* ~ = 2" (r) - ^'M = = "T * (T) ' 

{ftW.^W} = {O'M.^fT)} = (3.5) 

By going to Dirac brackets for the second-class constraints the Grassmann momenta are 
eliminated and we get (to simplify the notation we denote {., .}* with {., .}) 

{0«>(r),0«>>)} = -iSv, {el Q \r),ef\r)} = = 0, 

{^(t),0*(t)} = -iSv, {^(r),^(r)} = R*(r), 0*(r)} = 0. (3.6) 

The primary constraints are 

7r V(.)( T '^) ~ °' Tn(T,5) « 0, 7T n(a) (r,CT) « 0, 

3 M (a) (r, a) = 6(a)(6)(c) 3 e ( 6) r (r, a) 3 tt[ c) (t, a) « 0, 

7r r (r,a)^0. (3.7) 



By evaluating the canonical Hamiltonian by Legendre transformation (see Ref. [3]) and 
by asking that the primary constraints are constants of the motion under the r-evolution 
generated by it, we get the following secondary constraints 



16 Following Rcf . [33] , the classical regulated theory can be obtained by taking the mean value of the solutions 
of the final Hamilton equations with the Berezin-Marinov distribution function. For the U(l) group behind 
the charge rji = 8* 8i (topological two-levels for the sign of energy) the positive definite distribution function 
is pi — en + 8* 8i with a,i > 0, J pi ddi d8* = 1. Therefore the classical regulated value of a quantity A is 
< A >i2= / Api d8i dd\ p 2 d8 2 d6% where Eqs.(3.4) have to be used. 
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d r 7r r (r, a) + Qi Vi ^{°, Vi( T )) ~ 0, 



3„ 3_r 3_s 3n 

e e (o )( 6 )( c ) e (o) e (b) iZ rs(c) 



16ttG 



2ttG 

c 3 3 e 



3 Co(a)(6)(c)(d) 3 e( a )r ^(fi) ^^(cjs ^(d)] ( T > <?) + M (t, a) « 0, 



W( )(r , <?) = <9 r 3 7T(" a) - e( a )(6)( c ) 3 7T( c) + 3 e( a) .M r ] (r, <?) = 

= - 3 e[ a) (r, a)[ 3 e r - 3 M (6) ](r, a) » 0, 

3 r (r )( 7) = [ 3 7rf a) 9 r 3 e (o)i -9 a ( 3 e (o)r 3 7rf o) )](r,?)-^ r (r,a)w0. (3.8) 

In Eqs.(3.8) the following notation has been introduced for the matter terms (the mass 
density Ai(r, a) and the mass current density Ai r (r, a)) to conform with the treatment of 
the same matter in the non-inertial rest frames of Minkowski space-time given in Ref. [10] 



N 



M(t, a) 
Mi(r, a)c 

M r {r, a) 



* 3 (*, Vi(r)) Mi(r, a) c + 3 eT { ™\r, *), 



i=l 



Vi \Jm 2 iC 2 + 3 e r (a) (/c ir (r) - ^ Z) 3 e^ a) (kJj) - ^ A.^ (r, a), 

AT „ 

£ »* ("frOO - f A(r,a)) <5 3 (^(r)) - "eTiTV,*), 



2^ (4 3e(a)r 3e(a)s ^ ^ + y 3e w 3e w 3e w 3e w F ™ Fsv ) (T ' 

-i-F rs (r,a)7r s (r,a). (3.9) 

Let us remark that the mass current density M. r {r, a) does not depend upon the 4-metric. 

In Eqs.(3.8) T(t, a) « is the electro-magnetic Gauss law, 7i(r, a) « and H( a )(T~, o) ~ 
are the super-Hamiltonian and super-momentum constraints respectively, while 3 @(r, er) 
are the generators of 3-diffeomorphisms on E T . In the constraint 7Y(r, a) ~ we have 
3 e e (a)(6)( c ) 3e ( a ) 3 e(f,) ^rs(c) = 3 e 3 i? with 3 i? the scalar 3-curvature of the instantaneous 3- 
space S r and Q rs (a) the 3-field strength determined by the 3-spin connection ^uj r (a) given in 
Eq.(2.10). 

One can check that the constraints are all first class with the algebra (see Ref. [34] for the 
terms containing the Gauss law constraint) 
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{ 3 M (a) (r,a)/M (fe) (r,a)} = e (a )( b)(c ) 3 M (c) (r, 3)8* (ff, a ), 



{ 3 M (a) (r,a), 3 e r (r )C T)} = 3 M (a) (r,a) 



<9<5 3 (a,<r) 



fG r (r,a), a s (r )C ?')} = 3 6 r (r, a') ^- + 3 O s (r,a)^- 5 3 (a,a') 



da r 



St) 



3 u 3 m 3 d 3 n t? 
Vt n \ d/„\ t/u\ r„ 



'(a) c (a) ^(6) *(6) 



{t,B)T{t,B)5\B,<?), 



{H(T,a), 3 Q r (T,a')} = H(r,a') 
{H(r,a),H(r,a')} - 



^,d5 3 {a,a ) tt s (t, a) 3 , 

+ <Vs ^ T(T, 0-)t> (<7,CT ), 



<9a r ' ~'' s 3 e(r, a) 

3 e(a)( r ' ? ) ^(a)( r - s ) + 3e (a)( r ' ^W( r ' 
4o s e? a) [ 3 6 s + 3 u s{b) 3 M (b) ]j (r, a) + 



d5 3 (a,a) 



+ 



[ 3 s + 3 ^(6) 3 M (6) ]](r )( T')) 



V A d5 3 (a,a') 



da r 



(3.10) 



The energy-momentum tensor of the matter, T AB (r, a) 
put Mj c = rji Mi c) 



_2_ ,)S 



part-\-em 



S 4 g AB 



(r, a), is (we 



T rs (r,a) 



M(t, a) 



[ 3 e(l + n) 2 ](r, «?) ' 

3 e r {a) \{l + n) 3 el a) M s -n {a) M 



3 e(l + ny 



[ 3 e 



1 N 



3„r 3^s 



_7 7* e (a) e (fe) 

M 



((1 + n) 3 e^ (^ m (r) - 5* A m (r, a)) - n (a) M,) 



(1 + n) (K m (r) - ^ A n (r, a) j - n (b) M, j j (r, #(t)) + T ( 7 m) (r, a), 
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l r 1 



3 e 2 



+ E 



3 r 3 s 
e (a) e (fe) 



abcduv 



f ri (a) n (b) \ 



+ 4 



"-(a) n (6) n (c) 



3„ 3„ „u _-u i 
e( c ) u e ( d)t, 7T 7T I + 



+ E 

ablm 

+ E 



(1 +n) 4 

' e (a) H) "(a) M*>) + "(6) Ma) 



3 r 3 s 

e (a) e (6) 



abcdlmuv 



1 +n 



3 Z 3 m 

e (a ) e (6) 



1 



5ab 



EM 



3 m 

(e) e (e) 



3 m 3 1 n (c) ^(cQ \ 3 „ 3 „ 1 , , 



(3.11) 



With these conventions we have .M = 3 eT±± = 3 e(l + n) 2 T TT and .M r = — 3 eTj_ r = 

3 e (1 + n) Ea \a)r [ Es M^ ^ + "<«) ^ ■ 



1. T/ie Dirac Hamiltonian and the Weak ADM Energy 

In the rest-frame instant form of tetrad gravity [3] , after making the Legendre transforma- 
tion and after the addition of the DeWitt surface, the Dirac Hamiltonian has the following 
expression in terms of the weak ADM energy 17 in the canonical basis (2.7) 



H D = - E ADM + / d 6 a 



nH- n {a) H( a ) A T T (r, a) + A r (r) P r ADM + 



+ 



A n vr„ + A„ (o) n n(a) + \ na) n V(a) + fl (a) 3 M (a) + n vr 7 



(3.12) 



where A n (r, a), A n(a) (r, a), A^^r, <j), H( a )(j, a), A r (r), are the Dirac multipliers in front of 
the primary constraints and the mass density M. is defined in Eqs.(3.9). 

Since the scalar 3-curvature 3 e 3 R = 3 e e( )(6)( c ) 3 ej^ 3 e^ fi rs ( c ) of S T can be decomposed 
in the following way [3] 



3 e(r, a) 3 R(r, a) ^ S(r,a) +T(r,a), 

s(t, a) d = f [ 3 e M-) 3 4o ( 3 C 3 rL - 3 r« s 3 ri)] (r, a), 

rstw 

T(r, a) ^ £ 9 r ( 3 e V 3 9 ra (ft, 3 g vt - d t 3 g uv )) , (3.13) 



tuv 



17 Since we have Eadm = cP^ DM , [Hd\ = [E/c] = \mlt~ 1 ] and [t] = [I], Hamilton's equations are written 
as d T F = {F, Hjj}- To agree with the standard ADM equations we need a minus sign in front of the shift 
function. 
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it follows that the weak ADM energy of Eqs.(2.22) and the super-Hamiltonian constraint of 
Eqs.(3.8) can be written in the following form 



Eadm — c 



J d 3 a [M - 



s + 



+ 



2ttG 

c 3 3 e 



167T G 

)(6)(c)(d) "e(o)r "T(6) ~ e (c)s ^(d) 



E3^ 3 3 r 3 3 r 

Cj (a1IWrVrf1 e Cn.V Kw s 7T, 



abed 



M 



16nG 



(S + T) + 



2 n G ^ 3 3 3 r 3 3 
T^T" 2^ G ( o )(6)( c )(d) e (a)r 7T( 6) e (c)s 7r ( s d) 



c 3 3 e 



abed 



(r,a) « 0, 



1 A 

c 



Eadm + J d 3 a (nw)(r,£) 



J d 3 a 
2ti G 



[l+n)M (t,o)- 



16nG 



[l + n)S + nT (r,a) + 



d 3 a 



— (1 + n) 3 G (a)(b)(e)(d) 3 e (a )r \ b) 3 e (c)s 3 7rf d) 

abed 



(r,a). 



(3.14) 



B. The Hamilton Equations for the Particles 



Due to the presence of the Grassmann-valued signs of the energy, which are constants 
of motion like the Grassmann-valued electric charges ({^(t), H d }—0, {Qi(r), H D }—0), the 
particle Hamilton equations are defined in the following way 



ViVlir) = {vI(t),H d }, 
ViK ir (r) = ki(T) = {K ir (T),H D }. 



(3.15) 



The derived particle world-lines are (r) = z M (r, ffi(r)). Looking at the expression of the 
momenta K ir (r), i=l,..,N, we can define the following derived 4-momenta, corresponding to 
the ones of the standard manifestly covariant approach [11, 13], which satisfy the mass shell 
constraints for each particle also in the curved space-time M 4 in presence of interactions 
(however, as shown in Ref. [8], other definitions are possible for these derived quantities) 
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= (p\ 0) =rn^m 2 iC 2 + 3 e[ a) (r, ^(r)) K ir (r) 3 e s {a) (r, r^(r)) K is (r); 

pS a) = ^>^))^))> 

„M _ „(M) 4 

n A _ (m)4 _ 
Pi - Pi (E)^) - 

= (* - TT^M)-' = -(2^)(^ W ) + . < .,(n^,,^), 

Pi ^Pi =Pi ^)HPi =Pi ^sPi =m iC . (3.16) 



The Dirac Hamiltonian (3.12) generates the following Hamilton equations for the particles 
(see Eq.(3.3) for the first one) 



ViVK T ) = Vi e[a)( r >^( r )) 



(l+n)^ a) ( Kis (r)-^4) 



-n {a) 



Mi( r ))> 



d . .o_ ViQi \- U ( \ dA u(T,Vi(i~)) dA T (T,fJi(r)) 



Vi ~l~ K ir( T ) z 

CLT C 



+ 



Vi 3~ \r\i)i 



J~ i> ( T 



rrij c 



JJl + n) - 3 e (a )« 3 e (a)v (fi?(r) + 3 e u [b) n (6) ) (rjV(r) + 3 e v (c) n (c) 
1 9[ 3 e (a ) s 3 e( a)t ](T,^(r)) 



- (l+n) ^ (r ^ (r)) + 



MO")) 



Vl(j) + r4)^(b)](T,77i(T))J (#(t) + reJ c) n (c) ](r,77;(T)) 



+ 



g [ 3e (a) W(«)]OVftO0) 



e (6)s e (6)t ](r,77;(T)) ( ^(r) + [ e (c) n (c) ](r, Vi(r)) 



(3.17) 



where Ti r {j) denotes a set of relativistic forces, which in non-inertial frames of Minkowski 
space-time would be only inertial effects [10]. 

As a consequence, the second order form of the particle equations of motion implied by 
Eqs. (3.17) is 
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\ a) r\ a) s m iC {r,!(r)+ 3 el b) n {b) ) 



(r,77i(r)) = 



o ?7i Qi 



dA u (T,f)i(T)) dA r (T,f)i(T))\ ( dA T (T,rJi(r)) dA r (r, $(t)) 



+ 



<9t 



+ 



or 



rjirriiC 



dr 



\ 



o_ Tji Qi 3 3 



9A u (r,r/i(r)) <9A r (r, ^(r)) ' 



(r,^(r)) = 



+ 



dA T (T,ffi(T)) dA r (T,ffi(T)) 



drfi 



dr 



H 



drjf 



™* C 3e (a) 3e (a) 



1 ( d[\ e)s \ e)t \{T,m) ^( T ) + p e . Cf)n(/)](T ^ (T)) ) (^(r) + [V (9) n (9) ](r,^(r))) 



+ 



g [ 3e (e)^(e)](^^(^)) 



dr 



<9 6(e)r C( e ) u N 



<9^ 



M(r)) i)?(r) + [ 3 e ^n w ](r,^(r)) 



(3.18) 



Here Ti r {j) is the new form of the relativistic forces. 

If, as in Ref. [10], we define the non-inertial electric and magnetic fields in the form 
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def ( dK_ _ dA , 

1 (hi! <); 1 ' 



def 1 



(3.19) 



the homogeneous Maxwell equations, allowing the introduction of the electro-magnetic po- 
tentials, have the standard inertial form e ruv d u B v — 0, e ruv d u E v + ^ = 0. 

Then also Eqs.(3.18) take the following form 



d 



dr 

o f]i Qi 



3 e (a ) r 3 e (a)s rriiC (j)t( T ) + 3e (b) n w) 



\ 



\j (l+n) 2 -*e ic)u \ c)v (#(r) + 3 e^n (d) ) (#(r) + ^ n {e) ) J 
[E r + e ruv rji(r)B v ] (r, ffi(r)) + ^ ^ r (r). 



(3.20) 



C. The Hamilton Equations for the Electro-Magnetic Field 

The Hamilton-Dirac equations for the electro-magnetic field are 



d 



A T (r,a) = cn(r,a), 



dr 

— A r (r,a) = ^—A T + ^— 3 e {a)r s e {a)s ir s + 3 e s {a) n ia) F sr )(T,a), 
d 



— n r (r,a) = £ Vi Q^r) 5 3 (a, #(t)) + 



+ (^7 [(1 + n) 3 e 3 e\ a) 3 e\ a) 3 e^ % } F uv - n {a) { 3 e\ a) ir r - % } tt*)] ) 



(r, a). 
(3.21) 



1. Maxwell Equations 
Eqs.(3.21) imply 



n'(r, a) 



1 + n 



3,3 (a) 3,3 (a) (Frr ~ ? ' ' &\ a ) 77 '(a) ^vr) 



= - V- 4 ^(r, a) 4 g TA (r, a) 4 g sB (r, a) F AB (r, a). 



(3.22) 
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If we introduce the charge density p, the charge current density j r and the total charge 

Qtot = J2i Qi 011 S r 



1 N 

f (r,a) = 3 V 7/i Qi # (r ) <5 3 (a, ff t (r) ) , 

=>• Qtot = / d 3 a 3 e(r, a) p(r, a) = ^ Q 4 
^ i=i 

then the last of Eqs.(3.22) can be rewritten in form 
d 



da 1 



ir r (r,a) « - 3 e(r, a) p(r, s ), 



(97r r (r, a) o d 



dr 



da s 



(3.23) 



V- e (a) e (a) e (fe) e (6) - "(a) ( e (a) 7T - e (o) 7T J [T,a) + 



+ 3 e(r,a)T(r,a). 



(3.24) 



If we introduce the 4-current density s A (r, a) 



N 



s T (r,a) = 



s t, a = 



1 N 



(3.25) 



and we use (3.22), then Eqs.(3.24) can be rewritten as the standard manifestly covariant 
Maxwell equations in general relativity 

1 d 



, , . , , , V~ 4 9(r, a) 4 g AB (r, a) 4 g CD (r, ff) F BD (r, a) =-s c (r,a). (3.26) 
Eqs.(3.26) imply the following continuity equation due to the skew-symmetry of Fab 



d 



y/-*g(T,ff) da< 



= 0, 



or 



d 



V / 7(t^) p(t, a) 



1 t~ ~ ~~~r~ [ 3 e(r,a)T(r,a)] = 0, 



3 e(r, <r) <9<r r 



(3.27) 
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so that consistently we recover Q to t = 0. 

2. The Radiation Gauge for the Electro-Magnetic Field in Non-Inertial Frames in ADM Tetrad 
Gravity. 

In Ref.[10] there is a general discussion about the non-covariant decomposition of the 
vector potential A(t, cr u ) and its conjugate momentum 7r(r, a u ) (the electric field) into lon- 
gitudinal and transverse parts in absence of matter. Only with this decomposition we can 
define a Shanmugadhasan canonical transformation adapted to the two first class constraints 
generating electro-magnetic gauge transformations and identify the physical degrees of free- 
dom (Dirac observables) of the electro-magnetic field without sources. This method identifies 
the radiation gauge as the natural one from the point of view of constraint theory. When 
there are charged particles, this method allows to find the expression of the mutual Coulomb 
interaction among the charges in the admissible non-inertial frames of Minkowski space-time. 
As said in Appendix B of the first paper of the final version of Ref. [10] the reduction to the 
radiation gauge has to be done after fixing the 3+1 splitting of space-time and the associ- 
ated radar 4-coordinates, i.e. after fixing the gauge variables 9' 1 (t, a) and 3 K(t,B). In the 
next Section the radiation gauge will be used in the restriction of the Hamilton equations 
to suitable gauges. 

Here we extend the results of Refs.[10] to our class of space-times. 

If A = ^ r dl is the non-covariant flat Laplacian, associated to the asymptotic Minkowski 
metric and acting in the instantaneous non-Euclidean 3-space E T , its inverse defines the 
following non-covariant distribution 

1 11 

— S 3 (a, a) = c(a, a) = =, (3.28) 

with 5 3 (a, a ) being the delta function for S r . 

Then we can define the following non-covariant decomposition of the vector potential 
and its conjugate momentum (T(r,a) ~ is the Gauss law of Eqs.(3.8), Ar] em (T, a) = 
5 rs d r A s {r 1 <?), i] em describes a Coulomb cloud of longitudinal photons, see Ref. [35]) 

A r (r,a) = A ±r (r,a) - d r r] em (r,a), 

n r (r,a) = n r ± (r, a) + 5 rs d s f d 3 a' c(a, a') (T(r, a') - £ Q t Vi 5 3 (a' , jfc(r))) , 

Vem(r, a) = - J d 3 a' c(a, a) (V s d r A s (t, afj , {r? em (r, a), T(r, a)} = 5 3 (a, a), 
A ±r (r,a) = 5 ru Pr(ff)A s (r,a), n r ± (r,a) = ]T P r ± s (ff) tt s (t, a), (3.29) 



where we introduced the projector P[ s (cr) = 5 rs — 5 ru 5 s 



du d v 
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If we introduce the following new Coulomb-dressed momenta for the particles 

c or\ i 

=^ K ir (r) - —A r (T,ffi(T)) = k ir (r) - — A ±r (T,ffi(T)) 
c c 

we arrive at the following non-covariant Shanmugadhasan canonical transformation 



(3.30) 



A r (r, a) 


Vl(r) 




A±r(r, a) r] em (r,a) 


Vl(r) 






i — > 






7r r (r, a) 


/tj r (r) 




n r ± (T,a) r(r,a)^0 


ki r (r) 



{A ±r (r, a), 7ri(r, a)} = c5 ru P*>{3) 5 3 (a, a'), 
{V r i(r),^is(r)} = S r s S i:j . 



(3.31) 



The non-covariant radiation gauge is defined by adding the gauge fixing i] em (r, a) rs 0. 
As shown in Ref.[10], the r-constancy, dVem ^'°~ ) = {r] ern (r, a), H D } « 0, of this gauge fixing 
implies the secondary gauge fixing for the primary constraint 7r r (r, <?) ps 



d 



+ 



da' r 

1 + n(r, <?')) [ 3 e ((l)r 3 e (a)s ] (r, a) 
3 e(r, c?') 



A r (r,a)~- / d 3 a'c(a,a') —— [ 3 e? } n (o) ](r, a) F sr {r, a) 



dc{a,fjj{T) ) 
da' n 



(3.32) 



If we eliminate the electro-magnetic variables A T , n T , f] em , T by going to Dirac Brackets 
(still denoted {.,.} for simplicity), we remain with only the transverse fields A± r and n 7 ± 

(F rs = d r A ±s - d s A ±r ). 

Let us remark that in the radiation gauge the non-inertial magnetic field of Eqs.(3.19) is 
transverse: B r . = e ruv d u Aj_ v . But the non-inertial electric field E r = —F Tr = — d T Aj_ r + 
d r A T is not transverse: it has E± r = — d T A± r as a transverse component. Instead the 
transverse quantity is 7r^_ (it coincides with 5 rs E± s only in the inertial frames of Minkowski 
space-time), whose expression in terms of the electric and magnetic fields is 7i r ± (T,a) = 

rfi M«) H) ( E > - e ™ nU B v)} ^ °) ~ 6rS £i Qi Vi d s c(a, fjiir)). 

The electromagnetic part of the Hamiltonian (3.12) can be expressed in terms of the new 
canonical variables, since from Ref.[10] we have (n r = 3 e r ^ n^) 
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J d 3 a ( 3 - 



;i + n)T| e _ L m) +n r Tj e r m) ])(r, ( T) = 



n r f 3 eT| e r m) + W r )](r,a), (3.33) 



where the new energy-momentum tensor has the form 



e i T ' ^ J _L_L l T >°V - [ 7r ± 7r ± + ^ e (a) e (a) e (b) 6 (b) ? ru ? sv^ (T, a) , 



C(a)r C(a)s 



3„r 3_s 3_u 3_i> 



^t^T^t,*) = -F rs (r,a)7r s ± (r,a). 



(3.34) 



In Eq.(3.33) we rewrote the non-inertial Coulomb potential W(r) (a function of the 
particle 3-coordinates t)1(t)), found in Ref. [9], in the form 



W(r) = J d 3 a[(l+n)W {n) +n r W r ](r,a), 



W r (r,a) 



3 e ( a)r 3 e ( a) s ( _ r ;m ^ n ^ dc(a,ffi(r)) 



27r r ± -5 rm Q 



iVi 



2c 3 e \ x ^ ' da m 

dc(a,f)j(T))- 



6 ^ ^ — — J ( r ' a ) ' 



(3.35) 



In Minkowski space-time its limit to inertial frames is the standard Coulomb potential 



E 



Let us remark that now this potential depends not only on the particle 



i^j 4Tr\ffi(T)-ffj(T) 

positions, but also on the non-inertial electric and magnetic fields, on the 3-metric on £ r 
and on the lapse and shift functions. 

The energy-momentum tensor (3.11) can be written in the radiation gauge by means of 
the substitutions A r i— > A± r and 7r r i — > 7r^_ — r\i Qi ac (^( T )) ; Kir \— > k ir . 



3. The Dime Hamiltonian and the Hamilton Equations in the Radiation Gauge 

After the elimination of the variables i] em , T, A T , 7r T by going to Dirac brackets, the final 
form of the Dirac Hamiltonian (3.12) in the radiation gauge is 
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H D = - Eadm + I d a 



+ 



nH -n {a )H( a ) (r, a) + X r (r) P r ADM + 
A n vr n + A„ (o) vr n(a) + A^ {a) + /i (a) 3 M (a) (r,a), 



(3.36) 



with Eadm, H(r,a) and 7i( a )(r, <r) given by Eqs.(3.8) and (3.14) but with .M(r, a) and 
•M r (r, a) replaced with the following quantities 



M(t, a) 
Mi(r, B)c 



N 

Vi(r)) Mi(r, a) c + (^rfr ) (r, a) + W H (r, 5), 



i* ,k c* + 0^/3 ^ Q -2 Vm Vga (~^) -^A ±r ) (Mr) - ^ A ±s ) (r, a), 

y ars 

AT ~ 

M r (r, *) = J> (^r(r) - ^ A ±r (r, a)) 5 3 (<t, #(t)) - ( 3 eT| e r m) ) (r, a) - W r (r, a), 



+ 



(f>f^ )y j(r, a) = r^(r,a) [1 £ Q= V„ V« tt! ttJ 

= [r V3 E ^Ka V M (2 7T1 - J™ g Q iVi 9C{B d ^ T)) ) 



5 sn Qi m 



dc(a,ffj(T)) 



da r ° 



(t, a), 



s 



da* 



(3.37) 



Let us remark that also A / i r (r, a) does not depend upon the 4-metric like Ai r (T,a) of 
Eqs.(3.9). We have already given the various quantities in the York canonical basis. 

In the radiation gauge the Hamilton-Dirac equations (3.17) for the particles are replaced 
by the following equations 
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m 2 c 2 



:(r,^(r)) - 



(3.38) 



The quantity W(r) of Eq.(3.35) is a functional of the 3-coordinates fji(r) (and also of 
the 4- metric and the non-inertial electric and magnetic fields), which replaces the two- 
body Coulomb potential of the inertial frames of Minkowski space-time. In Ref.[10] it 

is shown that we can write rjiQi £ r (r, ^(r)) = -rjiQi dA±r{ gf i{T)) +r] i Q i ( dA ^ T & 

n dA ±r (T,ffi(r)) 8W(t) _ n F / -+ / \ \ dW(r) 
'h^i q t grf 'h J ~'±r\ l ■, 'fry 1 ) ) q^t ■ 

The first of Eqs.(3.38) can be inverted to get 



, 3 e (a )r 3 e {a)s h)f (r) + 3 e^ } n {b) j 
Viktor) = r] i m i c[—= = 

y (l + ^ - 3 e (c) « 3 e (c) „ (i)?(r) + 3 efo n (<0 ) (i)?(r) + 3 e^ e) n (e) ) 

+ ^A ±r (r,r/,(r)). 

c 



(3.39) 



In the radiation gauge the Hamilton equations for the transverse electro-magnetic fields 



A± r and n 1 , are 



d T A ±r (r,a) = {A ±r (r, a), H D } = 

= E 4. p T{S) [ < 1+ ">^»- 3 ^ - e ^ £ ft * ^P) + 



+ n (a) e v (a) F vu \(r,a), 
d T ir r ± (T,a) = {ir r ± (T,a),H D } = 



(1 + n) E s 3 eL M^) 



^ c 2 + E„, fe 3 e? fe) (M^) - f 3 efo (^(r) - f A ±t> ) 
"(o) (T,?7i(r))+ 
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+ 



(1 + n) E ( 3 e %) % } ( 3 4a) 3 ^) - s efo %)) d n F sv + 



svb 



_|_ f) r 3 p 3 P s 3 P V ( 3 P n 3 P m _ 3 P m 3 n \ p \ , 

-t- o n e ( e^) j _r s „y -t- 

-I- r) n \ 3 p 3 P s 3 p« ( 3 P n 3 P m _ 3 P m 3 p™ "| TP _|_ 



svb 



I ^ / 3 ^ n Pi *rr m I ^) 3 ^ n -rr- m ^) 3 ^ m ^" 

+ n (a) y e (a) o» n 7t_l + d n e (a) 7t_l - d n e (a) tt ± 

- E 3e w ^ - 5 « ^ n E * a ^0 r))) 



- E( 3e w^- 3e w^)E^^ 



g 2 c(a,r^(r))) 
d o l d o n 



- d n n [a) E ( 3 e? a) - ^ O E * Qi dC{ %^ T))) ] *) ) • ( 3 - 4 °) 



D. The Final Form of the Constraints in the York Canonical Basis. 

1. The Super- Momentum Constraints. 

If -D r ( a )(6) = o~abd r + 3 ^ r (a)(b) = o~abd r + £(o)(6)(c) 3 ^r(c) is the covariant derivative defined 
in terms of the spin connection given in Eq.(2.10), the super-momentum constraint (3.8), 
with M. given in Eq.(3.37), has the following expression [6] in the electro-magnetic radiation 
gauge 



w w = E 5r 37r w + E 3 ^(«)w 3?r w + E 3e w ^« 



d - E ^( a )( 6 ) + E 3e (" a ) = 



r6 



= E [ 6 * b ^ + E (3 ^ V ™ V « b - ^ <5a 1 V rb Ka] 0" 1 <9 U + 

rb u 

+ E [ 7 " a Q a V ra V u b ~ lab Qb Vrb Ka] 9 U R a + 

a 

+ \ [Qb Qa 1 Ka (0 r V ub - 8 U V rb ) ~ QaQ^ Kb (<9 r V ua - 8 U V ra )^ + 
+ ^ E Qb 1 Qw (Via Vvb — V u b Ka) Ktu K» j 
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+ 



V c / tun -VbV/ 7 

+ t^Q- 1 V vaM v »0. (3.41) 



See Eq.(Bl) of Appendix B for its explicit expression. 

The formal solution for 7rf(r, <r) = ^ ^ XLtab ^ ^ Q» Q^ 1 e<ab C(a)(b)) (t, a) given 
in Eq.(3.7) of Ref. [6] (see its Appendix D for the needed Green function) has only to be 
modified by adding the matter term in the last line of Eq.(3.41). 

In terms of the shear the super-momentum constraints take the following form (the func- 

(0) 

tional derivatives of the super-momentum constraint must be done before expressing 7r 4 in 
terms of the shear 0C( o )(&)|a^&) 



%) (t, B) d = f ( A- ( .) W %) + £ %) M v ) (r, a) = 

rb v 

rb u 

+ £ [^aaQaQ^ 1 V ra V ub ~ ~fab Q b V rb V ua U R n + 

a 

+ \ [Qo Qb" 1 Ka (<9 r Kb - <9 U Kb) -QaQb 1 Kb (d r V ua - d u V ra ) 

+ - £ Qa, 1 Qb 1 Qw (Kta Kb — Kb Ka) Kui ^ Kto) 



.3 



0- 1/3 Q^ 1 Kb £ 7a>n c - + 2 / 3 (Q~ b l Kb ^ - ^ E ^c" 1 K c <T( 6 )(c))] + 

c c 

+ <T 1/3 1 £ ^ (r, a) » 0. (3.42) 



2. The Super-Hamiltonian Constraint and the Weak ADM Energy 



The weak ADM energy of Eqs.(3.14) has the following form in the York canonical basis 
and in the electro-magnetic radiation gauge with M. given in Eq.(3.37) 
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E 



ADM 



+ 2 E 

abtwiuvj 



^abt ^abu Vwt ^iw Vvu Bj v 7T, 



QaQ^-QbQa 1 

3 e 3 R-T = 1/3 3 i? - 8 1/6 A 1/6 - T. 



(3.43) 



The functions iS(r, B) and T(t,B) are given in Eqs.(B8) and (B13), respectively 

In Eq.(3.43) we have also given the expression 18 of the scalar 3-curvature 3 R = 
3 R[(f), n , Ra] in terms of 3 -R = 3 R[9 n , Ra] and A, which are the 3-curvature of S T and the 
Laplace-Beltrami operator for the conformal 3-metric 3 g rs (detg rs = 1), respectively, and are 
given in Eqs.(B20) and (B12). This expression is needed to put [6] the super-Hamiltonian 
constraint in the form of the Lichnerowitz equation for the conformal factor = 1//6 , which 
has the following expression in the York canonical basis ( 3 g rs is the inverse of the 3-metric 
with unit determinant) 



H(t,B) 



167T G 

2tt G 



(r, B) [8 A 1 / 6 - 1 / 6 3 R ] (r, B) + M (r, B) + 



+ 2 E 

abtwiuvj 



T/ R T/ R ^) 
^abt ^abu * wt ^iw * vu ^jv ^Tj 



(t,<t) « 0, 



QaQ^ 1 - QbQa 1 ^ 

A = d r ( 3 g rs d s ) = 3 g rs 3 V r 3 V S = 9 r (E Qa" 2 K-a Ka 9.) = 

a 

= E Qa 2 [Ka V sa <9 r d s + (2 V ra V sa E 75a 9r R~b ~ 9 r (Ka Ka)) 9 a ] • (3.44) 



In terms of the shear the super-Hamiltonian constraint and the weak ADM energy take 
the following forms (again these forms cannot be used for the functional derivatives of these 
quantities) 



H(t,B) 



1nG r 
8tv 2 G 2 



■\r,B) (^(A- 1 - 3 R)^+ 2 -^^M + 

l2ir 2 G 2 ~ 



1 ~, 



o /I Lt ^ — \ « 1 72 \ — v 



ab,a^=b 



2 nf)(T,B)^0, 



18 



See Ref.[26] and Eq.(B4) of Ref. [6]. 
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E 



ADM — C 



ab,a^b 

The super-momentum and super-Hamiltonian constraints are coupled equations for <fi and 

(0) ~ 

7Tj- (or (f> cr ( a ){b)\a^b) ■ in the form given in Eqs.(3.42) and (3.45) the constraints have to be 
solved to get the conformal factor and the off-diagonal terms of the shear 0"( o )(&) \a^b- 

3. The Dirac Hamiltonian in the Radiation Gauge 

The Dirac Hamiltonian (3.12) takes the following form when written in the York canonical 
basis and restricted to the electro-magnetic radiation gauge 



H D — - Eadm + 



n 



U - n {a ) W( fl ) (t, B) + A r (r) P r ADM + 



+ / dV 



An vr n + A n(o) 7r n(a) + A^j 71-^^ + A a{a ) 7r^ (r, <r), 



(3.46) 



where A Q (r, a) are new Dirac multipliers replacing the /X( a )(r, (?) appearing in Eq.(3.12) 



E. The Asymptotic ADM Poincare' Algebra, the Rest-Frame Conditions and the 
Center of Mass in the York Canonical basis 

While the weak ADM energy is given in Eq.(3.43), the other weak ADM generators of 
Eqs.(2.22) have the following form 



pr 

r ADM 



-2 Jd 3 a{- 0- 2 / 3 Qf [VraV va ( ^rjr^ + ^TbaTaa ~ (1/2)^) 3 v R a ) - 



-^Tbad v (V ra V va ) 



a,v,b 

n 5 + 



+^ 1/3 J2Qa 2 [VraV va ^r 1 d v 4> + J2^ad v R a ) - ^d v (V ra V va )] 7T^ + 

a,v a 

+4>~ 2/3 Q« W [^«Kd ( + Yl 7aa d v Ra) ~ \v rd d v V va + 

a,d,v a 

+\ J2QfQlV rc V uc V vd (d v V ua - d u V va )] J2 n e n V : tB ni-i ~ 

-^" 2/3 E VraV sa M s } » 0, 
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J ADM = Jd 3 a{a r [- 0" 2 / 3 Qf [VsaV va ( Irba^dJ + ^(TbaTaa ~ (1/2)^) ^) 



n- h + 



a,v a 
+0 _2/3 E [VsaVvd^-'dJ + ^2laa O v R n ) ~ ^V sd d v V va + 

+^ E Q7 2 Ql v scV U cV vd (d v V ua - d u V va )j E 



1 7 



2/3 ' 



t - Q d Q-' - QaQa 1 

r J2Qa 2 VsaV ua M u ~ 

a,u 

-v s [ - V 2/3 E ( + E(r 6a 7,a - (i/2)^ 5 ) ^ik) - 



a,u a 

4 Y,Qc 2 QlV rc V uc V vd (d v V ua - d u V va )} E n ^^ tjg nn-i " 

A c.u twi ^ d ^a ~ ^ a ^d 



J2Qa 2 VraV ua M u ] }, 
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J 



ADM 



(-' „ , . 2ttG' i _ 1 



16nG 



S-M- 



^ 1 (-3(^) 2 + 2^n? + 



+2 E 



VJD T/ D (^) (^) 



)] 



+ 



abtwiuvj ^QaQb 1 _ QbQ a 1 

^0 V3 E {<5a 2 [KaKa( - + ^ 7aa + ^(KaKa) 



+ 



a,v 

1 7 



+ 



+ (Qa 4 - E ^c" 2 )^(KaKa)] } } « 0. (3.47) 



The rest-frame conditions P\ DM ~ and the conditions J^dm ~ eliminate the internal 
3-center of mass of the 3-universe E T . The gauge fixings J^dm ~ imply A r (r) = in 
i^u of Eq.(3.36), because their time preservation implies <9 r J^dm — {Jadm-> ^d} = Padm ~ 
\ t {t) Eadm ~ —\-(j)Eadm ~ 19 . In this way we identify the Fokker-Pryce center of 
inertia of the 3-universe as the origin of the 3-coordinates in the instantaneous 3-spaces S T . 



F. The Dirac Hamiltonian in the Schwinger Time Gauges and in the Electro- 
Magnetic Radiation Gauge 

In what follows we restrict ourselves to the Schwinger time gauges ar a \(T,a) ~ 0, 
(/?( a )(r, (?) ~ 0, whose r-preservation implies A v . At, a) = X a , At, a) = in Eq.(3.46). The 

following results are obtained after the elimination of the variables «(„), 7r^ and </?( a ), 7r (/ j (o) 
with Dirac brackets. 

By using Eqs.(3.44) and (Bl) for the super-Hamiltonian and super-momentum constraints 
with Ai and Ai r given in Eqs.(3.37), the Dirac Hamiltonian (3.34) in Schwinger time gauges 
and in the electro-magnetic radiation gauge is 



19 The asymptotic Poincare' charges are assumed gauge invariant, i.e. {^AW'^t 7 "^)} 

{J T ADM^{a){ T , #)} W 0. 
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H D = - Eadm + I d a 



nH-n {a) H (a) (r,a) + 



/ 



+ / d 3 a 



A «7r +A (a) vr (a) 



= J d 3 a 
2txG 



[l + n)M 



16nG 



d 3 a\(l + n)S + nT 



J, <?) + 



J d 3 a [(1 + n) r 1 ( - 3 (0 tt^) 2 + 2 ^ n- 2 + 



a^b 



+ 2 E 



E e * £a(ra Kit o (0) ->\ 
-2 2^ 5 ™ 2^ s j« n ) J J ( r ' 

abtwuv \Q a Q b — QbQ^ 1 ] i 3 



/j3 - / -»\ 1-1/3/ ->\ ( tabtQf, 1 V r bV w t B iw (q\ 

b^a rtwi ^ b ^ <^ a ^ 

+ E(E 



tabt 



Qfe 1 &r (Yrb Vwt B iw ) + 



+ 



— TZl a n n -l ^ (7ca - 7cb) <9r -Rc Kb -B iw + 

+ E E n ^-^n-i ^ V ™ ~ V ™ V ^ * V " b V ^ B ™) ^ + 



+ Qa 1 E V ™ + E 75a<9 r n 5 ) + 

r b 

+ 1 E 9r Vra ~ Vra d r R b + Yl Tbb Vua Vrb 9r Vub ) U ~ b + 

rb ub 

+ Q- 1 J2 V ™ jCi r)(T,a) + 

r 

+ J d 3 a jA n 7r n + Af?( a ) 7Tfl( g )j (t, a), 



(3.48) 



where we have used the last expression in Eq.(Bl) for the super-momentum constraints. 

Hp depends not only upon the tidal variables R a , Ha and the matter, but also upon the 
gauge variables n, n^, 9 n , n^, which play the role of inertial potentials. 

Let us remark that in the Hamiltonian (3.48), and in particular in the weak ADM 
energy (3.43), the kinetic term / d 3 a 0(1 + n) 7r|] (t, 3) = -J^ I d ^ a H 1 + 

n) (^K^ ( r ) 3) connected to the momentum gauge variable 7t^(t, <x) = ^ G 3 K(t, a), de- 
termining the instantaneous 3-space, is definite negative in every gauge. Therefore it plays 
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the role of a dark energy and it vanishes only in the CMC gauges 3 K(r, a) ~ 0. Instead, the 
kinetic term connected to the momenta n\ (r, a) (or by the off-diagonal terms of the shear 
&(a){b)\a^b(T, <?)), determined by the super-momentum constraints, is always positive definite 
due to Eq.(3.45). 
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IV. THE EQUATIONS OF MOTION IN SCHWINGER TIME GAUGES. 



In this Section we shall write the Hamilton equations for the gravitational and matter 
variables, generated by the Dirac Hamiltonian (3.48), in an arbitrary Schwinger time gauge 
and in the York canonical basis. 

Let us remark that in the evaluation of the Hamilton equations the constraints can be 
used only after having done the Poisson brackets with Hp. When, like in Section V, some 
gauge fixings will be added, we shall restrict the Hamilton equations of this Section to the 
chosen gauge and we will not restrict the Dirac Hamiltonian to the gauge and evaluate the 
new Hamilton equations with the new Hamiltonian (even if the two approaches should be 
equivalent). 

By comparison in Appendix A there are the standard ADM equations of canonical metric 
gravity and a discussion of the contracted Bianchi identities. The use of the York canonical 
basis allows to disentangle the contracted Bianchi identities from the equations for the 
remaining gauge variables (the inertial effects), for the tidal variables and for the matter. 

See Appendix B for the explicit form in the York canonical basis of many terms appearing 
in the Hamilton equations. 



A. The Contracted Bianchi Identities 

The variables 0(r, a) and nf\r, a) are the quantities determined by the (non-hyperbolic) 
partial differential equations corresponding to the super-Hamiltonian and super-momentum 
constraints, respectively. The Hamilton equations for them are the contracted Bianchi iden- 
tities ensuring the r-preservation of the constraint sub-manifold: they hold independently 
from the form of the solution of such partial differential equations. 

For the Hamilton-Dirac equations obtained by using the Dirac Hamiltonian (3.48) in 
the York canonical basis, in the electro-magnetic radiation gauge and with the Fokker-Pryce 
observer as origin of the 3-coordinates are 



d T 4>(T,a) = {4>(T,a),H D } 



5H 



D 



5 71 At, a) 



\2txG 



+ n {a) 



(1 + n) H + 2/3 Qa 1 (dr n (a) V ra + 

ra 

) + d r V ra ) (r, a) 



2 

Vra [t;^ 1 9 r <p-^2 Tba 9 r R b) + d r K 



where Eq.(B2) was used. 



(4.1) 



For Tif^ the Hamilton-Dirac equations obtained by using the Dirac Hamiltonian (3.48) in 
the York canonical basis are 
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d T nf\r,a) = {nf\r, a), H D } 



(»), 



SH. 



D 



se^r, a) 

■5M{t,B x ) c 3 6S(T,ai, 



--/^[P + »)^(7^-«C5^)- 



16ttG 

47TG 



//-, . \ 7-1 V"^ tabt tabu n i 77 j ^ [Kut -Siw Bj v ] \ . 

( (1 + n)0 *» ) M + 



+ 



f ^3 - / 



r,<7 



(4.2) 



where Eq.(B21), (Bll), (B19) and (B5) have to be used. With Eq.(2.18) we can replace vrf ) 

With <7( a )( 6 )| a ^ 6 . 

B. The Equations of Motion for the Gauge Variables 

The equation of motion for the lapse and shift functions identify their r-derivatives with 
the arbitrary Dirac multipliers 



d T n(r,a) = A n (r,a), 
d T n {a) (T,a) = Afi (o) (r,a). 



(4.3) 



For 7r^ the Hamilton-Dirac equations obtained by using the Dirac Hamiltonian (3.48) in 
the York canonical basis give the following form of the Raychaudhuri equation 



drni(T,a) = {7T7(r,a),iJ D } 



5H 



D 



l,- 5(r , s) 



S(j)(r,a) 6 r 5(j)(r, a) 



c 



+ 



16ttG 
2tt G 



n{r, a x ) 



ST^ax) 



+ 



c 

a^b 



(l + n)(T,?) 

tabt tabu Vwt K 



3^ + 20- 2 (^nf + 



aftiumu j^Qa Qb 1 ~ Qb Q a *] 

+ (<T 1/3 n (o) g^ 1 9 r tt^) (r, a), 



(4.4) 
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where Eq.(B3) has been used and Eqs.(B22), (BIO) and (B17) are needed. With Eq.(2.18) 
we can replace itf^ 1 with <J( a ){b)\a^b- 

For 9 % the Hamilton-Dirac equations obtained by using the Dirac Hamiltonian (3.48) in 
the York canonical basis are 

dr = {9\r,a),H D } * H ° 



57rf\r, a) 

8tY G , \ 7 i ^abt ^abu Vwt Ku ^iw -Bjv ^ j 



abtwuvj (Q a Q f} 1 — Qt> Q a 1 j 

- 0- 1/3 (r,a)E [n (o) (EE[^ _1 ^ + 

a bj^a rtw 

„ 1 CaW Qb 1 K-6 Kot -Bju) . 

+ £ £ brfQa ^^L_ (y rc \/ ua _ \/ m \/ uc ) Q r y ub y wf B iw ) — 



e bct Q a 1 QbQc 1 

I 

bu c^b 

tabt Qb 1 K& Kit -Bj 



where Eq.(B4) has been used. 

Then one has to add the gauge fixing constraints (satisfying the orbit condition) to the 
super-hamiltonian and super- momentum constraints 



x(t,<t) « 0, 



(clock synchronization convention or 

determination of and of the instantaneous 3 — space) 



Xr(r,a) « 0, 

(determination of 9 n and of the 3 — coordinates) . (4.6) 



Their preservation in time, by using the Dirac Hamiltonian (3.48), generates the equations 
for the lapse and shift functions consistently with the clock synchronization convention and 
with the 3-coordinates (the partial r-derivatives act on the possible explicit r-dependence 
of the gauge fixings) 
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or 

(determination ofn) , 



drXr(r,a) = dX f ,a) +{Xr(r,a),H D } = 0, 
a r 

(determination ofn,( a )), (4.7) 

and as a consequence the Dirac multipliers are determined: A n (r, a) = d T n(r, a), 
Ag ( a)( r > #) = d T n {a) (T, a). 

For gauge fixings of the type 1 (t,(j) ~ (numerical function) 1 and n^(r,a) m 
numerical function, Eqs. (4.7) are just Eqs.(4.4) and (4.5), respectively. 

C. The Equations of Motion for the Tidal Variables and the Matter 

Let us now consider the Hamilton-Dirac equations of motion implied by the Dirac Hamil- 
tonian (3.48) for the tidal degrees of freedom and for the particles. 



1. The Hamilton Equations for the Tidal Variables 

For the tidal variables R a we get the following kinematical Hamilton equations 

d T R a (r,a) = {R a ( r , a), H D } = = 
8nG ~ x 



(i + n)n, 

c 6 



- <T V3 E Qa 1 (»(«) [Taa V ra (- ^ d r + £ Tba d r - 

ra b 

- V ra d r R a + ^2 ^~ ab Vsa Vrb dr Vsb ] ~ 7aa dr n ( a ) Vr °) ( r > ( 4 - 8 ) 

sb 

where Eq.(B6) has been used. 

Eq.(4.8) can be inverted to get the momenta Tig, in terms of the velocities d T Ra 
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n s (r,a) = 



d T R- a + 



8ttG l + n(r,a) 
+ <r 1/3 (n {a) [ Taa V ra (i 0" 1 9 r $ + T ba d r Ri) 

ra I 

- V ra d r Ra + ^2 ^ b Vsa Vrb ^ Vsb ~ ^ aa ^ r n ( a ) Vr > 



sb 



The dynamical Hamilton equations for U a are 



(4.9) 



d T U- a (r,a) = {U- a (r,a),H D } = 



SH. 



D 



SR a (T, a) 



a) 16tt G 5 Ra(r, a) 



16nG 
8tiG 



n(r, <7i) 



.JT(r,5i) 



5i?a(r, a) J 



c 



^ ^ (laa lab) 



^abt tabu Vwt V vu B% w Bj v (Qb Q a ~\~ Qa Qb , 



abtwuvij 



{QbQa'-QaQb 1 



+ J dViti (a) (r, ai) 



_ <JW( a )(r,ai) 



5i?a(r, ct) 



*f *f )(r,*) + 



(4.10) 



where Eqs.(B23), (B9), (B15) and (B7) have to be used. With Eq.(2.18) we can replace vrf ) 

With <7( a )(6)|a#>- 

If we evaluate the r-derivative of Eq.(4.9) and we equate it to Eq.(4.10), we get the 
following second order equation for R n 
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d 2 T R- a (T,a) = [ E Qa'Vran^YsiTaaTba - 5-ab)d r d T R- b + 

ra I 

+ 4>~ 1/3 E [(> V ra (\ ^ ^ + E ^ ^ R c) - V ra R~a + 

ra c 

+ E lab V sa V rb d r Kb) n(a) - 7aa K« d r n (a) ] E 7 5a <9 T i? 5 - 

sb 6 

- <T 1 [flr i? s + I V 1/3 E 1 ( [7*. Ka {~\ ^ r 4> + E Tba ^ i%) 

ra ft 

- Ka <9 r + E 7ab Ka Kb 9 r Kb] W( a ) - 'Jaa V ra 8 r W( ) ) ] <9 T - 

sb 

- ^ _4/3 E 7aa 1 Ka "(a) #r <9 T 0- 



- <T 1/3 E ^ [ " ("(a) - ^ (5 ^ ^ + 

ra 

+ E ^ 5r ^)] + 7sa dr ™(°)) + E 756 

b sb 



d V sa V rb d r V sb 



e T e l + 



+ [d T R- a + 0" 1 / 3 E ^a 1 ( [l-aa V ra (- 0" 1 «9 r + E 7ab <9 r - 

ra 5 

- V ra d r R a + E 7ab Ka Kb <9 r Kb] "(a) ~ 7aa Ka <9 r W( a )) ] 



<9 T n 
1 + n 



- 4>' 1/3 E <5a 1 ( [laa Vra (\ V ' r + E Tba 9 r R~b) ~ Ka r R~ a + 

ra b 

+ E 7sb ^ a Vrb d r V sb d T n {a) - 7„ a V ra d r d T n (a) ) (r, a) + 



+ i (1 + „)) (r, 9) / [(1 + „) (T , *) m| + „(„ , |gg 



a^b 

( E (^aa-7afe) 
abtwuvij 



^abt tabu Kut Ku -Bjui -Sju (Qb Qa Q« Qb / 



*fM*>)(r,*). 



(QbQ^-QaQ^ 1 )' 



(4.11) 
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This equation depends: 

a) on the r-derivatives <9 r 0(r, <r) and d T 6 l (T,a), given in Eqs.(4.1) and (4.5), respectively; 

b) on the r-derivatives of the lapse and shift functions, namely on the arbitrary Dirac 
multipliers appearing in Eqs.(4.3). 

2. The Hamilton Equations for the Particles 

By using Eqs.(3.15), the Hamilton equations (3.38) for the particles in the electro- 
magnetic radiation gauge take the following form in the York canonical basis (Eq.(3.35) 
and (3.37) are used for W) 



ViVKr) = M(r),H D } 



dMjr, a) 

9 K/j r 



= J d 3 cr (jl + n(T,a)j 

a v lr 

0" 4 (1 + U) Z as Qa 2 Vra V sa (k ts {r) - & A ±s ) 

= Vi [ i = 
Jmfc 2 + 0" 4 J2 CUV Q~ 2 V uc V vb (k iu {r)-^A ± ^ (k iv {r) - & A ±v } 

- <T 2 Yl Qa 1V rafl(a)Y T ^i( T ))^ 



rii-^Kirir) = {Kir(r),H D } 



= - J d 3 cr (Jl +n(r,a)] 



dMjr, a) 
drjl 



ST- i ^ L-2 Y^n-iT/ 1/ dM v (r,a) 



)- 



d 



/^[(l + n)^ + 



+ f a x;Q. _i fi(.)E 7 « 



dr}] 
dW v i 



dril 



r, <?) + 

(1 + n) Q- 2 V sa V va (k iv {T) - f A ±v } 



^ 0- 2 g: 2 K a n (a) )f ^](r,r/,(r)) + 



+ ^F ir (r,r/i(r)), 
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V cmn 



Qi 



^im(r) A± m i j />,,, 

C 



Kin I T 



c J 



On 



(««(r)-^i s ) 



+ 



mf c 2 



<P~ A Ecmn Qc 2 VmcV nb {k im {r) - ^ A ±m j («i„(r) - ^M ±n ) 

+ 0- 4 (l + n)£ 2 Vua (V sa (\ ^ O r + £ Tn , <9 r ^) " 9 r Ka) 

(« itt (r) - % A ±tt ) (k is (r) - & A ±s ) 
(^m? c 2 + 0" 4 £ cnm Qc" 2 Kn C Kb (K im (r) - f Auj (Mr) - f A ±n ) ) 
- 0- 2 ]T g^ 1 n (a) (V sa (1 r 1 9 r + £ 75. 5 r R- b ) - d r V sa ) 

as I 

(k is (r)-^A ±s ) 

yQ C 2 ~ -4 g-2 y mc V nb (k im ( T ) ~^A ±m ) [k m {r) ~%A ±n ) 

By inverting the first of Eqs.(4.12) we get the following form of Eq.(3.39) 



[r,Vi(r)). 
(4.12) 



Ki r ( T ) 



A ±r (T,f) i (T)) + 



- 2/3 E ^cKcKc (i)r(r) + 0" 1/3 £ Q, 1 V«,n w ) (i)?(r) + 0" 1/3 £ Q; 1 V ve n (e) )] ^ 

mic d e 

M(r)). (4.13) 



If we put Eq.(4.13) into the second of Eqs.(4.12), we get the second order equation for 
t}1{t) (corresponding to Eq.(3.17)) 
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^=(-^W + — Vt(r) ^ + Vi F„ ) (r, *(r)), 



0tf 



W(„)(r,5) = 



W(r) = / dV [(1 + n) W (n) + 0" 2 ^ g^ 1 n (a) ^ Ka W„ 

a u 

4 [*- E « K» (* -1 - ^ E Q. n. ^P) 

ars \ i / 



(7, <T), 



W r (r,a) = --F rs (r,?)rVQ !?|i 



Fir C 



1 + ra - 



J] Qc Kc V vc (^(r) + 0- 1/3 £ Q d 1 



-1/2 



a 6 

+ 4 Ql V ™ {{\ ^ d r j> + J2 Taa O r R a ) V va + d r V va ) 

auv a 

(vUt) + r 2 £ Qb 1 ^n ( b)) (#(t) + 0- 2 £ g," 1 K c n (c) ) 



(4.14) 



5. TTte Hamilton Equations for the Transverse Electro- Magnetic Field 



Eqs.(3.40) for the transverse electro-magnetic fields A ±r (r, a) and 7t^_(t, a) in the radia- 
tion gauge take the following form in the York canonical basis 
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d T A ±r (r,a) = {A ±r (r, a), H D } 



iVi 



dc {a,rji{r)) 



= S m PF(?) [<T 1/3 (1 + n) Ql Vua V a -^rJ^Q 

nuva m i 

d T n r ± (T,a) = {ir r ± (T,a),H D } = 

wma i 

r ^ 2l \l + n)Q~ 2 V ma Y. s V sa k ls {r) 



- r 1/3 Q- 1 V ma n (a) ](T,ff l (r)) + 



+ 



+ 



(1 + n) (^" 1/3 2 ^ 2 Vsa V vb (V na V mb - V nb V ma ) d n F sv + 

svbn 

+ 0~ 1/3 Q- 2 Qf [d n (V sa V vb (V na V mb - V nb V ma )) - 

- Vsa V vb (V na V mb - V nb V ma ) 0- 1 d n + 2 Y (Tba + lib) d n R b )] F sv ^j + 

b 

+ 0" 1 / 3 9 n n Qa 2 Qb 2 Vsa V vb (V na V mb - V nb V ma ) F sv + 



svbn 



+ 



+ 



-1/3 



9 n V na ~ Ka ( ■= 4> 1 9 n + ^ Tba 9 n Rb 



7T 



- \d n V ma ~ V ma (- 1 <9 n + ^ Tba ^n R~b 
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V n a V na 



1 d n 4> + Y J rMdnR i )\s"" - 



- [d n V ma - V ma (i 0" 1 9 n + £ 75a d n R b )] 5 nt ) £ Vi Qi dC{ff d ^ T))) - 

t i 

- 0^ £ dn U (a) ^ (Ka <T* - Kna O £ V i Qi <?) ) • 



(4.15) 
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4- The Constraints to be added to the Hamilton Equations 

To the previous Hamilton equations we have to add the super-Hamiltonian constraint 
(3.44) [or (3.45)] for the determination of 0(r, a) (the Lichnerowicz equation) and the super- 
momentum constraints (3.41) [or (3.42)] for the determination of the momenta nf\r, a) [or 
0(t, <?) <7( a )(b) \a^b( T , <?) of Eq.(2.18)] or their explicit expression given in Ref.[6]. 

Also the six rest-frame conditions contained in Eqs.(3.47) have to be added. 
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V. GAUGES IN ADM CANONICAL GRAVITY IN THE YORK CANONICAL 
BASIS 



Let us consider some gauge fixings for the secondary super-hamiltonian and super- 
momentum first class constraints determining the gauge variables vr^(r, a) and 9 n (r, a) inside 
the family of Schwinger time gauges. 

As shown in Eqs.(4.7),their preservation in time by using Hp will generate four secondary 
gauge fixings determining the lapse and shift functions consistently with the chosen definition 
of instantaneous 3-space (clock synchronization convention) and of 3-coordinates in it. 



A. ADM 4-Coordinate Gauges 

These are CMC (constant mean curvature) gauges with 3 K(r, a) ps and with the 
following gauge fixings for the 3-coordinates 

n^r,a) « 0, 

Xr{r, a) = ["^2d s 3 g rs - ^d r ^ 3 9ss (t, a) = 

s s 
a s 

+ 2j2(l^ 1 ds4> + Ylr b ad s R- b )V ra V sa - 

s b 
~ \ dr^ + Y,rbadrRl))]{T^) « 0. (5.1) 



These are 3 equations of the form Ylas ^ia( T > °) ®s T( rs ) a (9 n (r, a)) + 
Eas F 2as(r,B)T {rs)a (e n (r,a)) + F 3 (r,a) = with T (rs)a = V ra {6 n ) V sa {9 n ) = T (sr)a 
Once T( rs ) a (r, a) has been found, one has to find the compatible Euler angles n (r, a) of the 
rotation matrix V(9 n ). 

A similar gauge with 3 K(r,a) ks and with Xr(T,&) ~ replaced by 3 g rs {T,(j) ~ [(1 + 
|<^) 4 <5 rs + h^J\{r, a) (ti^J is a transverse-traceless tensor) is used in Ref.[36]. 

B. The 3-Harmonic Gauges 

These are gauges in which only the 3-coordinates are fixed but not the instantaneous 
3-space 
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7r^(r, a) ~ not specif ied, 

s s a 

Y,Qa 2 (d s (V ra V sa ) + 



+ {\r 1 dj-2 1 £ Tba d s Ri) Vra V sa )] (r, 3) » 0. (5.2) 

6 

Again we have a linear partial differential equation for T( rs ) a = V ra (9 n ) V sa (9 n ) = T( sr ) a . 
C. The 4-Harmonic Gauges 

The 4-harmonic gauge, see Eq.(6.4) of Ref. [6], is defined by the four gauge fixings 
X A (r,ff) = Eb 9b (v^V B )M » 0. By using Eqs.(2.6), (2.10), (4.1) (the kine- 

matical Hamilton equation for 0) and Eqs.(4.3) (the r-derivatives of the lapse and shift 
functions are Dirac multipliers) we get the following expression for x r ( r i c?) ~ in the York 
canonical basis 



X' 



12nG 



d T ((1 + n) 3 eVj +J> ((1 + n) 3 eVj 



d T — — d s 6n 1 = e d T — d s — — 

N ^ N J I 1+n ^ 1 



63 e( a) n (a) 



+ n 



4 



7T7(r,a) = 3 ir(r,a) 



(l + n(r,5)) 



^——2 (An + V 1/3 E ^ ^« ^ n ) ( r > *)> 



0, 



(5.3) 



The other three gauge fixings are 
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X r = d T ((1 + n) 3 e V ) +J2d s ((l + n) 3 e V') 
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]- e ^9 s ((l+n)0 6 ^ 3 e[ a) 3 e 
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1 +n 

.2 ~ 



£ «' ( 



"(a) + 
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d T e i + 



V ra (- 0" 1 <9 r - ]T 75a 5r - V ra — =-1 ) 

o ^-r' 1 + nJ / 

Qt 1 (VraV sa ((1+n) 2 5 

(ib + "-(a) ™(b) ) d s n - 



[l + n)< 



sab 



1 + n) [ n (a) d s n (b) + n (b) <9 S n (a) 



+ 



+ (1 + n) ((1 + n) 5 ab - n( ) n(6) J [<9 S (Ka Kb) + 

+ Ka Kb (^ 0- 1 <9 S - ]T ( 75 a + 75b) d s R- b )] ) ) (r, a) » 0. 



(5.4) 



The kinematical Hamilton equations (4.1), (4.5) and (4.8) are needed to get this final 
form. 

These unconventional Hamiltonian constraints (x T ~ does not define a CMC gauge 
3 K(r,a) = const.) are four coupled equations for 7r^ and 6 l in terms of 0, i?a, ILi, n, 
\ n = d T n, n( a ), A f j( a ) = d T n^ a ). 

The stability of these gauge fixings requires to impose d T Xaij, <?) ~ and <9 T x r (r, <?) ~ 0. 
In this way we get four equations for the determination of n and n^ a y But these are not 
equations of the "elliptic" type like with ordinary gauge fixings. They are coupled equations 
depending upon n, d r n, d T n, <9 2 n and ni a \ , d r nt a \ , d T nt a \ , S 2 . ri( a ) , namely hyperbolic equa- 
tions. As a consequence there is a problem of initial conditions not only for R a but also for 
the lapse and shift functions of the harmonic gauge. Each possible set of initial values should 
correspond to a different completely fixed harmonic gauge, since once we have a solution for 
n and ?!(„) the corresponding Dirac multipliers are determined by taking their r-derivative. 

Instead of reading these constraints as gauge fixings determining 7r^(r, a), Eq.(5.3), and 
9 n (r,a), Eq.(5.4), let us solve them for the Dirac multipliers A n and A S(a) . Then it is 
more natural the following interpretation: the family of 4-harmonic gauges is determined by 
all the Hamiltonian gauge fixings of the type (4.6) for vr^(r, a) and 6 n {r,a), with induced 
secondary gauge fixings (4.7) for the lapse and shift functions, such that the Dirac multipliers 
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A n (r, a) = d T n(r, a), \ n(a) (r, a) = d T n^(r, a), have the form implied by Eqs. (5.3) and (5.4). 



As a consequence, let us remark that in the family of harmonic gauges there is no natural 
way to visualize the instantaneous 3-spaces (they cannot be Euclidean due to the equiva- 
lence principle) and to check which properties have an inertial origin due to the freedom 
in the clock synchronization convention (the gauge freedom in the York time). Most of the 
applications (for instance the IAU conventions for the solar system [37]) are based on post- 
Newtonian expansions inside a fixed Euclidean 3-space of an inertial frame in Minkowski 
space-time containing Newton gravity as the zero order! 

D. The Synchronous gauges 

The synchronous gauges are defined by the gauge fixings 



implying A g Tr (T, a) ps (no gravito-magnetism) and 4 g Tr (r, <r) xs e ( 1 + n(r, <r) 



These unconventional gauges have a residual gauge freedom in any fixation of the 3- 
coordinates implying A s . a )(r, a) ~ 0. If we add the "comoving" condition n{r,a) ps 0, we 
are also restricting the freedom in the fixation of the York time. 

These gauges are used in cosmology in presence of Killing vectors implying that the solu- 
tion of Einstein's equations are homogeneous and isotropic. However their use in absence of 
Killing symmetries is questionable, since gravito-magnetism cannot be eliminated in general. 



n (a)(r, B) ps 0, 



(5.5) 
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VI. THE 3-ORTHOGONAL SCHWINGER TIME GAUGES. 



Let us now consider the most natural family of Schwinger time gauges in the framework 
of the York canonical basis, i.e. the 3-orthogonal gauges 

6\t,g) « 0, (6.1) 

in which the 3-metric is diagonal in each point of the instantaneous 3-space S r and V ra (8 n ) = 
5 ar . The residual gauge freedom is the selection of the instantaneous 3-spaces by fixing the 
gauge variable 7r^(r, a). 

Till now we have not make any hypothesis on the type of angles 6 1 (t,(t) to be used 
to parametrize the matrices V ra (9 l (T, <?)) and its derivatives V^y a = 9Vr d a jji ^ |g'=0; ^(i)jw — 

9Bj ^gf - |e.=o needed in the 3-orthogonal gauge. 

It is convenient to parametrize the group manifold of the Lie group 0(3) in terms of 
canonical coordinates of first kind (see Appendix A of Ref.[38] and references therein), 
because then we have: 

A) 6 n are first kind coordinates (not the standard Euler angles); 

B) V(6 n ) = e^*» " IT \ where (Tj) a6 = e iab (T; are the generators of the 0(3) Lie algebra); 

C) V(i)o6 = Ciab, #(i)a6 = \ e iab with B = A^ 1 , A^ ab = \ ti a b- 

In Appendix C there the restriction at the 3-orthogonal gauges of the quantities evaluated 
in Appendix B. 

In the 3-orthogonal gauges Eqs.(2.18) become 

c 3 



7T {d) \ 

71 i \9*=0 



c 



a 

3 



8ttG 



^ £iab<t><T(a,){b)\c&bQaQb \ 



ab 



7 87T G \- e a fej 7T, 

0O-(a)(6)|a^ft = "g- 



c 3 Y QbQ^-QaQ?' 

AttG 

7 3 



3 K rs \ et=0 = 2/3 Qr <5s 0-(r)( S ) |r^s - ^^^rsQr ( 2 1 7ar Ilg - 7T^ . (6.2) 



A. The Weak ADM Energy 

The weak ADM energy of Eq.(3.45) becomes 
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EadmIo^O — c J d 3 (T M\oi=0 — -^Q n Q ^le'=o + 

+ ^0-i(-3(^) 2 + 2£n! + 2£ 



W Jf>) 



Cabi Cabj 7Tj 7T 



= c J d 3 a [M\ 



abij ^Q a Q b 1 ~ Q& Q a 1 



6» l =0 



167rG 



•5 U<=n + 



E 2 67rG I 21/ -\ 



(6.3) 



where Eq.(6.2) has been used and Eq. (C8) has to be used. 

The mass density appearing in Eq.(6.3) has the following expression implied by Eqs.(3.37) 



M(r,a)\ ei=0 = Vi(r)) V i Jm* + ^ £ Q~ 2 (k ia {r) - ^ A ±a ) * (r, a) + 

i y a 



T]_ ± \t, a) 



= <T 4/3 (r, a) [1 & 5 ™ *~ < ^ + £ E 2 <?6 2 ^ ^] 

rsa aft 

W (B ,(r, a) = -1 [0-/3 £ g2 ^ ^ ( 2 ^ _ £ ^ £ „ 



dc{a,f}j{T)) 



(6.4) 



B. The Super-Hamiltonian and Super-Momentum Constraints 



From Eqs.(3.44) (or (3.45)) and (3.42) the constraints determining and 7rf^ are 
(Eqs.(C12), (C20) and (6.3) are also needed) 
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12 7T 2 G 2 



a) (6) 



ab,a^b 



i 12 7r 2 )(r,a)^0, 



(6.5) 



^abi Q a 



E (^ca - leb) d b R z vrf } + 
t, a) = 



b^a i lQ b Qa 1 - QaQ b j c 

+ Q- 1 (0 6 «9 a ^ + E (r 6a d a u- b - d a R- b n- 6 ) + A4 (o) 

= -^^ 2/3 m(E^ _1 h ff ww + 

+ <9 fe + E (75a - 7&b) 9 6 #5) <7( a )( 6) ] - 

b 

~ Qa 1 [^d a ^ + Y^(TbadaTlb - d a Rb^b) + M (a) ])(r,a) «0. 

6 

(6.6) 



To get the second form of the super-momentum constraints we used Eq.(6.2), which implies 

d r 7ff V=0 = gf^ Ea^fe € iab Qa Qb~ l 9 r <7 (a)(b) + (^T 1 <9 r + J^b (iba ~ Tbb) d r R b ) C(a)(6)] • 

For the solution of Eq.(6.6) the results of Ref.[6] are needed. 



C. The Contracted Bianchi Identities 



From Eqs. (4.1) and (4.2) the Hamilton equations for the unknowns in the constraints 
are (the first is the Raychaudhuri equation) 



66 



\2ltG 



co (l + n)0^ + 2 / 3 ^Q- 1 (9 a n (a ) + 

a 



(6.7) 



d T nf\r, a)\ 6 i 



- J dVi 

.3 



c 



167rG 
8ttG 



/-, w - x /SM(t, <j\) . 
5 T(t, o\ 



c 3 5«S(T,CTi). 



V <50*(r,a) lw ' =u 16ttG <S0*(t,£) W= V 



se^r, a) 



a+b / T/ D \ (0) (0) 

n , u-i T e abteabj{V(i)kt + n k vr • ■ 



+ J d am^iT,^) §0i ^^ \o*=o, 



4 



9 T a (a)(b) \a^b,9 l =0 



87TG 



^abi 



QbQa 1 ~QaQ b - 1J ' 

5 .M (r, <7i 



r, a 



dVi [(l + n)(r,5i) 



c 3 6S(t,(Ti 



<50*(t, a) 



16nG 59 i 



r , Q\) 1 \ 
r,(j) / 



16ttG 



n(r, ai) 



((1+^)0 E 



C(m)(n) 



^ ^ £ mni (V(j)H + -B(j)fe t ) ffecd <5c Qrf 1 C(c)(d) ) ( r , &) 

5H( a )(r, cti) 



+ 



J3 



+ / d <7i7l( a )(T,<7i) 



50*(r, a) 



(6.8) 



Eqs.(C21), (Cll), (C19) and (C5) are needed. 



D. The Shift Functions 

By using Eq.(4.5) for the r-preservation of the gauge fixings (6.1), we get the following 
equations for the shift functions (the second expression uses Eqs.(6.2)) 
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a 

^2 1/3 (1 +n) o-( o) (6) | a#6 ] (t,<t). 
This is the final form of the equations for the shift functions. 



(6.9) 



(6.10) 



E. The Instantaneous 3-Space and the Lapse Functions 



If we recall from Eq.(2.9) the definition = j^-q K, we can restrict ourselves to the 
family of gauges 



F(r,a) » 0, 



^(r,a)^F(r )C T), 



(6.11) 



where F(t, a) is a numerical function independent from the canonical variables, which de- 
scribes the relativistic inertial effects associated with the freedom in the clock synchroniza- 
tion (they do not exist in Newton theory in Galilei space-times). 

For F(r, a) = const, we have the CMC gauges. 

By using Eqs.(4.4) and (6.2) restricted to the 3-orthogonal gauges, the preservation in r 
of the gauge fixing constraint (6.11) gives the following equation for the determination of 
the lapse function 
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where Eqs. (C22), (CIO) and (C17) have to be used. 

Note that this equation depends also upon the shift functions. 

Eqs. (6. 10) and (6.12) are 4 coupled equations for the lapse and shift functions. 



(6.12) 



F. The Equations of Motion 

1. The Tidal Variables 

The second order equations of motion(4.11) for the tidal variables R a and Eq.(4.9) for 
the momentum U a become respectively 
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where Eqs. (4.3), (6.7), (C9), (C15), (C23) and (C7) have to be used. 

The last term in the equation for d 2 . R a can be written as — ^(1 + n) E 



2. The Particles 



The Hamilton equations (4.12) for the particles become 
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The inversion of the first of Eqs.(6.11) gives 
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5. The Transverse Electro- Magnetic Field 

Finally the Hamilton equations (4.15) for the transverse electro-magnetic fields in the 
radiation gauge become 
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VII. CONCLUSIONS 



In this paper we gave the Hamilton equations of ADM tetrad gravity coupled to charged 
scalar particles and to the electro-magnetic field in the York canonical basis first in an 
arbitrary Schwinger time gauge and then in a sub-family of 3-orthogonal gauges. The 
electro-magnetic field has been specialized to the non-covariant radiation gauge, in which 
we can eliminate the electro-magnetic gauge variables and visualize the non-inertial Coulomb 
potential among the particles. 

We gave some refinements of the York canonical basis of Ref. [6], connected with the 
congruence of the Eulerian observers associated to the 3+1 splitting of space-time. 

We emphasized the role of the inertial gauge variable 3 K(t, a), the York time, connected 
to the freedom in the choice of the non-dynamical part of the clock synchronization conven- 
tion defining the instantaneous 3-spaces: in the York canonical basis it is the only gauge 
variable described by a momentum (a reflex of the Lorentz signature of space-time) and it 
gives rise to a negative kinetic term in the weak ADM energy vanishing only in the gauges 
3 K{r,a)=0. 

Moreover, in connection with the weak ADM Poincare' charges, we showed that in these 
asymptotically Minkowskian space-times one can introduce an interpretation of the isolated 
system gravitational field plus matter like in the inertial and non-inertial rest frames of 
Minkowski space-time [8, 10]. The 3-universe contained in an instantaneous 3-space can be 
described as a decoupled non-covariant non-observable pseudo-particle carrying a pole-dipole 
structure, whose mass and spin are those identifying the configuration of the "gravitational 
field plus matter" isolated system present in the 3-universe. 

In the next paper [19], starting from the results obtained in the family of non-harmonic 
3-orthogonal Schwinger gauges with an arbitrary numerical function describing the inertial 
York time, 3 K(t, a) ps F(t, a), we will define a linearization of ADM canonical tetrad grav- 
ity plus matter, to obtain a formulation of Hamiltonian post-Minkowskian gravity (without 
post-Newtonian expansions) with non-flat Riemannian 3-spaces and asymptotic Minkowski 
background: i.e. with a decomposition of the 4-metric tending to the asymptotic Minkowski 
metric at spatial infinity, a Qab 4 Vab- We can write a Qab = a Vab + 4 h,AB, with the small 
perturbation 4 /iab vanishing at spatial infinity, but this decomposition has no intrinsic mean- 
ing in the bulk differently from what happens in harmonic gauges with a fixed background 
(in them one has Euclidean 3-spaces with a violation of the equivalence principle). We will 
show that a consequence of this approach is the possibility of describing part (or maybe 
all) dark matter as a relativistic inertial effect determined by the gauge variable 3 K(t, a r ) 
(not existing in Newtonian gravity, where the Euclidean 3-space is an absolute notion): the 
rotation curves of galaxies would then experimentally determine a preferred choice of the 
instantaneous 3-spaces. 

Finally in a future paper we will replace the matter with a perfect fluid (for instance 
dust [25]) and we will try to see whether the York canonical basis can help in developing 
the back-reaction [39] approach to dark energy. 
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APPENDIX A: THE COMPARISON WITH THE STANDARD ADM EQUA- 
TIONS OF MOTION, THE RAYCHAUDHURI EQUATION AND THE CON- 
TRACTED BIANCHI IDENTITIES. 



In absence of matter, after the 3+1 splitting of space-time and the introduction of the 
radar 4-coordinates, the standard ADM equations of motion of canonical metric gravity are 
(see Eqs.(4.11) of Ref.[2]) 

d T 3 g rs = n r \ s + n s \ r - 2 (1 + n) 3 K rs , 
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-U 3 



(Al) 



Since from Eqs.(2.6) and (2.14) we have 7 = det 3 g rs = ( 3 e) 2 = 12 = 2 and 3 K = 
7T7 {6 is the expansion of the non-geodetic Eulerian observers), we can deduce 
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The second of Eqs. (A2) is an equation for the time evolution d T ir^ of the gauge variable 7r^, 
named Raychaudhuri equation. For geodetic congruences of time-like curves it is connected 
with the geodesic deviation equation and is a fundamental ingredient for the singularity 
theorems, see Ref. [40], where it is shown that in general (also in Minkowski space-time) in 
a congruence of "time-like geodesies" caustics will develop if convergence (9 < 0) occurs 
anywhere. In certain space-times there will be real space-time singularities if other global 
assumptions hold [40]. As shown in Ref. [41] the singularity theorems imply that the gauge 
fixing identifying the instantaneous 3-spaces (clock synchronization convention) have to be 
divided into inequivalent classes. Those who satisfy the strong energy conditions and have 
non-negative 3-curvature imply a singularity (geodesic incompleteness) in the past if the 
spatial average of the expansion does not vanish. In our class of space-times, including the 
Christodoulou-Klainermann ones [42] in absence of matter, the assumed boundary conditions 
should eliminate the singularities (at least for finite intervals of time; but in presence of 
matter the problem is open). Moreover the congruence of Eulerian observers, relevant for 
the York canonical basis, is not geodesic. 

Instead the first of Eqs.(A2) is a contracted Bianchi identity. Let us clarify this point. 
Due to the Bianchi identities 4 G^ '.„ = 0, the 10 Einstein's equations G^ u = — 
\^9p-v^R= -^-qT^v imply the equations of motion T^ U ]U =0 for matter and 4 contracted 
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Bianchi identities, which say that 4 Einstein's equations are not independent from the others 
and their time-derivatives. Therefore there are only 6 Einstein's equations functionally inde- 
pendent. But, since 4 Einstein's equations are the super-Hamiltonian and super- momentum 
constraints (restrictions on the Cauchy data), it turns out that only 2 combinations of the 
10 Einstein equations are dynamical, i.e. depending on the accelerations (they can be put 
in normal form; see the Hamilton-Dirac equations for R a , Il s in Section III). The four con- 
tracted Bianchi identities, A G^ V = 0, imply [40] that, if the restrictions of Cauchy data 
are satisfied initially and the spatial equations = are satisfied everywhere, then the 
secondary constraints are satisfied also at later times. Behind these identities there is the 
gauge invariance of the Einstein-Hilbert action under 4-diffeomorphisms. However to get 
a canonical formulation we need a 3+1 splitting of space-time and the replacement of the 
Einstein-Hilbert action with the ADM action 

At the Hamiltonian level, while the first 6 equations (Al) are nothing else that the def- 
inition of the extrinsic curvature of the instantaneous 3-space S T , only two of the other 6 
equations (Al) are independent due to the ADM version of the 4 contracted Bianchi iden- 
tities, now implied by the gauge invariances of the ADM action which are no more the 
4-diffeomorphisms due to the presence of the 3+1 splitting of space-time. At the Hamilto- 
nian level these gauge invariances are generated by the 8 (14 in tetrad gravity 20 ) first class 
constraints of ADM canonical metric gravity defining a pre-symplectic sub-manifold of the 
phase space symplectic manifold. However at each fixed value of the radar time r the gauge 
group is not a Lie group 21 and its group manifold, as a transformation group on the symplec- 
tic manifold of ADM canonical metric gravity, does not have 8 functional gauge parameters 
but only 4, the 4 Dirac multipliers 22 in front of the 4 primary constrains (the vanishing 
of the lapse and shift momenta). The secondary constraints (the super-Hamiltonian and 
super-momentum constraints) have as gauge parameters the gauge variables correspond- 
ing to the lapse and shift functions, i.e. 4 elements of a Darboux basis of the symplectic 
manifold. As a consequence, it is not clear how to define an abstract group manifold for 
this type of symplectic gauge group 23 . Moreover, the group manifold is r-dependent: we 
have (a-priori) different group manifolds for the gauge group on each instantaneous 3-space 
S T . Therefore Eqs.(Al) contain an information on the r-dependence of the effective group 
manifold of the gauge transformation group of canonical gravity describing its modification 



In canonical tetrad gravity there are 6 extra primary constraints with the associated Dirac multipliers in 
the Dirac Hamiltonian and no extra contracted Bianchi identity. 

It contains spatial diffeomorphisms acting on S r and the gauge transformations, replacing the time diffco- 
morphisms, generated by the super-Hamiltonian constraint, which modify the shape of the instantaneous 
3-space Er with deformations along the normal to S T in every point, proportional to the conjugate gauge 
variable 7rr = 3 K, namely to the value of the York time in that point. 

They correspond to a special choice of the arbitrary generalized velocities existing at the Lagrangian level 
due to the gauge invariances. The kinematical set of the Hamilton-Dirac equations identify them as the 
r-derivatives of the lapse and shift functions. 

From what has been said the natural gauge parameters of the symplectic gauge group are the lapse 
and shift gauge functions and their r-derivatives, so that there are only 4 arbitrary functions replacing 
the space-time 4-diffeomorphisms (whose abstract group manifold, depending on 4 arbitrary space-time 
functions as gauge parameters, is not under mathematical control in the large), the gauge group of the 
Einstein-Hilbert action. 
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from an instantaneous 3-space £ r to a modified one S T +dr- 



Let us now look at the Hamiltonian interpretation of the ADM equations (Al) in the 
York canonical basis of tetrad gravity. Since in the canonical York basis there is an iden- 
tification both of the gauge variables (a( a ), V?(a)j n, fi( a ), 9 r , n^) and of the variables (0, 

(0) 

7Tr ) determined by the super-Hamiltonian and super- momentum constraints, the content 
of Eqs.(Al) rewritten in the York canonical basis will be: 

i) 4 physical (hyperbolic) equations for the tidal variables R a , ILj, see Eqs.(4.8) and (4.10). 

ii) 4 equations for the variables 0, see Eqs.(4.1) and (4.2), which hold independently 
from the form of the solution of the super-Hamiltonian and super-momentum constraints 
determining them. They are the Hamiltonian version of the 4 contracted Bianchi iden- 
tities implying the time-preservation of these secondary constraints. Let us remark that 
the super-Hamiltonian constraint (3.44) and the super- momentum constraints (3.41) are 4 
coupled elliptic partial differential equations for the unknowns an d with coefficients 
depending upon the gauge variables 6 l and (not upon the lapse and shift functions) and 
upon the tidal variables R a , n^. Therefore their general solution (depending also on some 
arbitrary functions) will be of the form F[R a , H a , @\ Kg], 7r] ~ Fi[R a , U a , 9\ -Kg]. 

iii) 4 equations (4.4) and (4.5) for the gauge variables 6 r , n^: therefore, after the addition 
of gauge fixings (4.6) these equations allow to evaluate the terms d T 9 r , d T ng appearing in 
Eqs.(4.7) and generate the equations for the determination of the lapse and shift functions. 

By comparison let us consider electro-magnetism (with its Abelian gauge group) in the 
Shanmugadhasan canonical basis (3.31). In the rest-frame instant form of Ref.[8] the canon- 
ical variables and the constraint are A T , n T ~ 0, r] em = —j^d-A, n Vem = T = d ■ if 0, Aj_, 

7r_i_. The Dirac Hamiltonian is H D = f d 3 a nj_ + B 2 — A T n Vem + A 7r r j (r, a) and the gauge 

variable A T plays the role of the lapse and shift functions (\ = d T A T ). The analogue of the 
equations (4.4), (4.5), for the gauge variables is d T r] em = — A T : this equation determines A T 
once the gauge fixing for i] em has been given (the analogue of Eqs.(4.7)): rj — f[..,A T ] m 0. 
The analogue of the Hamiltonian contracted Bianchi identities (4.1), (4.2), is d T n Vem =0 24 . 

In conclusion in the York canonical basis we have the following interpretation of Eqs.(Al): 

1) The first six equations (Al) involve the 3-metric 3 g rs = 2//33 g rs with det 3 g rs = 1 
and 3 g rs = J2 a Q\ V ra {6 1 ) V sa (9 l ). While 7 = det 3 g rs = 4> 2 is determined by the super- 
Hamiltonian constraint, 3 g rs = <\>~ 2 I 3 3 g rs depends upon the three gauge variables 9 r (the 
freedom in the 3-coordinates on E T ) and on the tidal variables Ra through Q a = 7sa-Ro_ 
Therefore the first six equations (Al) imply Eq.(4.1) for and 



24 The final equations for 4> an d in tetrad gravity are so complicated because the gauge algebra of 
metric gravity is non- Abelian and not a Lie algebra. 
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a a 

o_ 2 r 12ttG 
~ ~3 I &~ 



2 
3 

+ 0" 2/3 



1 + n) 7T7 



n ia) 



+ 



(0 2/3 £ V ra {0 1 ) n (a) ) s + (0 2/3 ^ g a Ka(^) n (a) ) - 2 (1 + n) 3 K r 



(A3) 



The five equations (A3) are a mixture of the two Hamilton-Dirac equations (4.8) for R a 
plus the three equations (4.5) for d T 9 l : in the York canonical basis it is possible to separate 
the two Eqs.(4.8) from the three Eqs.(4.5). 

2) The second group of six equations (Al) involve the extrinsic curvature tensor given 
in Eq.(2.10). These six equations split into the second one of Eqs.(A2), i.e. Eq.(4.4), for 
the trace of the extrinsic curvature plus five equations for the traceless extrinsic curvature 
3 K rs = 3 K rs — | 3 g rs 3 K, which, in the York canonical basis may be decoupled to give the 

two Eqs.(4.10) for H a and the three equations (4.2) for Ttf\ 
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APPENDIX B: CALCULATIONS. 

In this Appendix there are the calculations needed for the form in the York canonical 
basis of the Hamilton equations in arbitrary Schwinger time gauges, generated by the Dirac 
Hamiltonian (3.48), of Section IV. 

1. The Functional Derivatives of the Super-Momentum Constraints. 

a. The Super- Momentum Constraints 

For the evaluation of the functional derivatives it is more convenient to use the following 
form of the super- momentum constraints obtained from Eq.(3.41) 

■Zj I -\ .4-2/ ->\ ( €aht Qb 1 Vrb Vwt B iw iq\ 

n {a) (r,a) = (r, a) I 2^2^ n n-i-n n-i drn * + 

b+a rtwi ^ b ^« ^ a 

+ V ( V r—— t iQZ 1 d r (V rb V wt B iw ) + 

Tttt, b^a QbQa ~ Qa Qb ^ 

+ n n-i n n-i E ^ lsa ~ 7a ) dr Re Vrb Vtw Biw + 

Qb Q a Qa Qb c 

+ E E Q n Q n-i ^ V ™ ~ V ™ V ^ * V <* V * B ™) ^ + 

Qc Qb Qb Q 



bu c=£b 



+ Qa' E V ~ {f 9 r H + E Tba * 11*) + 

r b 

+ Qa 1 E { Tba 9r Vra ~ Vra dr R ~ b + E V " a 5r ^ ub ) 115 + 
rfc ub 

+ Qa l Y. V ™ M r){T,t) »0, (Bl) 

r 

This form of the super-momentum constraints will be used to evaluate their functional 
derivatives with respect to the arguments = 6 , n^, 9 l , irf\ R a , H a - 

b. The Functional Derivative with Respect to 7ri 



f j3 - / - x 5H( a ){T,ai) 





= "[^E^ 1 (?rfi {a) V ra + 
ra 

+ n (a) [K a (40- 1 ( 9 r 0-^75a«9r^)+-9rV;a])](r, ( ?). (B2) 
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c. The Functional Derivative with Respect to (j> = 6 
Since we have _ - / - = I </>~ 5 (t, a) S J ~ , we get 

<JW(o)(r,5i) 



rfVi n (o )(r,5i) 



5 0(r, a) 



6 [0 3 ]T n (a) Q- 1 K-a d r ttJ (r, a). (B3) 



d. TTte Functional Derivative with Respect to irf 1 



By using the results 
5, 



Q^ 1 



and 



Qb Qa 1 -Qa Q, 



— — 



Q Qa 1 QbQc 1 

r Qc Q b —Qb Qc 



QbQa -QaQu 



rz?5 o- 1 ^5 75b - {lla ~ lib) 



Qb Qa 1 +Qa Q b 1 
Qb Qa Qa Qu 



d r Rb 



QcQZ 1 +QbQ7 1 



we get 



f J3 / - X ^^(a)(r, CTi) 

rfff iL n (-)( r ' ^ T (9)7 ^ = 

a 

/ \ ^ L ,-1 n , , o n 1 e abtQh VrbV w t B 

■hi 



QiQa 1 -QaQb 1 



+ 



b=£a rtw 

bu c+b Q c Qb -QbQ c 

tabt Qb Vrb Kut 



QcQb —QbQ 

d r^) 



b^a tv 



QbQ- a l -Q a Q b 



[r,a). 



(B4) 



e. TTje Functional Derivative with Respect to l 

By using the result = g>(ff, E r ftr 



E r ^fgF^5 3 (^<?i),weget 
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a y ' 

-(f)~ 2 {t, a) Y d r n (a) (T,a) 



b sb 



+ E(E 

jtw by^a 

+ EE 



tabt Qb 1 9 V r b V w t Bj w 



QbQ- a l -QaQ~ b 



+ 



tbct Q a 1 QbQ 1 



r\ t r 

r (Kc V.a - Ka V.c) ^ K,t Tff ] (r, 5) + 



QcQb 1 — Qb Q c 



+ 2 (t,^) Yl n («)( r '^) 



+ ( 75(1 (2 0" 1 dr^ + Y, Tea 9r Rc) ^b ~ 9 r R~ b II 5 ) + M r 



~df T» [Vsa Vrb (d r U b - (2 0" 1 d r + ]T 7 ca «9 r R 5 ) II 5 ) + 



sbb 



+ d r (v m v rb ) n 5 ] + J2 iu dr v,„ n s ) - 



- E E 

rstwjb c^b 

+ E(E 



tbet Qa Qb Q c n/ T/ T/ T/ n dVsb T/ R o (0) , 
^ ( VVc ^sa - Vra ^sc j "T^- V wt B jw r 1T j + 



rst wj b c^b ^cQft -QbQ 
tabt Qh 1 



+ EE 

sb c^b 



, b*a Q^-QaQb 

ebctQa'QbQ- 1 -0\: 



90* 



QcQ^-q^q? 



d6 i 



(^r'drt + Y (r b a-Tbb + 



75c 



— <9 r 



QcQ; 1 -QbQ; 1 
(V rc V sa V^. a V^ c ) Bj w 

d (V rc V sa — V ra V sc ) Vwt Bj w 



(V rc V sa -V ra V sc )V wt B 



JW 



TV 



(r,a). 
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/. The Functional Derivative with Respect to Tl a 



f J* - / .^W(,)(t,5i) 

Jda^n {a) {r^) = 

ra b 
V ra d r R a + ^2 ^ab Vsa Vrb 9 r Kb] ~ laa 9 r W( ) V ro ) (t, (?) . 



g. T/ie Functional Derivative with Respect to Ra 

By using the results 

= n 7 , a g«(r,a)5 3 (a,a 1 ), 



<5fls(r,?) Q b Qa 1 -QaQ b 
<5 3 (<7,<7i), 

5 3 (a,<ri), 

<5 Q n 1 Qc 1 

5 3 (CT,<Ti), 
9, 



i QbQa 1 -QaQ b 1 



(QbQa^Qa O 2 



Q 



r _ r _ T ( r> ^ 1 ) = 



6 Vc 



(QbQa -QaQ^ ) 2 (Qb Qa -Qa Q b ) 

we get 



lab {laa lab) Q t Q~ 1 -Q a Q~ x 

7aa / I7aa 7afej Q^~q^I 
laa ~ {lab ~ lac) Q cQ ^-Q bQ -i 



= 0~ 2 (r, ct) <9 r n(o) (r, a) Q" 1 \/ ra II S - 

ra 



a 

1 ^ («9r Ka + Ka II a - (2 «9 r + ^ 7ba 9 r R b ) U- a 

r b 

- laa [ ^ (tSo #r K a + 7 56 V sa V rb d r V sb ^j U b + 

b sb 

+ v ra (<p 6 d r 7ct + Y(rb a d r n- b -d r R- b n- b )+M 
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Ta 



EE 



'"" Q7 ' - i V ri V„ t B l „d,.xf' + 



b^a rtwi 



' QbQa'-QaQb 



e abt Q a 



+ EE( 2 E 

i rtw by^a 
+ J2{laaTbb-labTba)drR- b 



-1 



(QbQ-'-QaQb 1 ) 2 

Vrb V w t B{ w 



2 (laa ~ lab) <f> 1 <9 r + 



laa ^ ^ 



^abt Qb l 



-EE 



e bct Qb Qc 1 



— 8 r (V r b V wt B iw ) 



QcQb 1 - QbQc 1 



laa ~ (lab ~ Ta 



2 Qa Qb 



b c+b 

( V ™ V «* - V ra Kc) 9 r V sb V wt B iw j 7ff ) (t, B) 



QcQk 1 - QbQc 1 ! 



(B7) 



2. The Function S and its Functional Derivatives. 



By using Eqs. (Bl) and (B2) of Ref.[4] we have the following expression for the function 
S(r,a) defined in Eq.(3.13) 

C d lf 3„ 3^ 3„s / 3 r u 3 r v 3 r u 3 r v \ _ 

° — e 2^ (o) («) v ™ sn _ rs u v ~~ 



= 2 E 

a 

( E [ 2 ( Vra d vln<p + 5 rv E Ka d u In <pj - 

rv u 

- 2j2v rb V ub V va Q 2 a Q b - 2 dJn ( p + 

bu 

+ E Tbb Vrb { 5ab 9v R ~ b + Vvb E V " a du Rb ) ~~ 

bb u 

- J2 Tba Vrb V ub V va Ql Qf d u R- b + 

bbu 

+ g E Vrb Vua [p u Vvb + dv Vub ) + 

+ \ Yl Vrb Vwb Q c ^ { 6ac [ d " Vwc ~ dw V ™] + 



bcw 



+ v vc y y u< 



f) v — 8 V 



) 



X 
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|2 (v va d r lri(p + 5 rv Y V sa d s In j - 

s 

- 2 J2 V vd V sd Vva Ql Qf d s ln<P + 

ds 

+ Y 7Ed Vvd { Sad dr R * + Vrd S ^ sa ^ 

cd s 

- Y ^ca K d V sd V ra Ql Qf d s R- c + 

cds 

+ \j2 V ^ V ^{ d ^ + d rVsd) + 
ds 

+ \j2 V ^ V ^eQd 2 (^ 



+ 



+ Ke J2 V ta 



d t V se - d s V te 



+ 6 2 Y Qf r 1 d r \d s (V ra V sa ) + 2 V ra V sa (0" 1 O s - Tea d s Re) 



= 2 Y Yl ( 2 ^ d r <t> [ d * ( V ™ V sa) + 2 V ra Vsa (2 0" 1 d s - ^ 7 c"a <9 S R 5 )] + 

a rs c 

+ Ka V sa Y( 2 75a 7ca ~ #r #5 <9 S - 



6c 



+ 



2 Y [rta Va d s V ra + Y Tbb Vra V vb V sb 8 V V 8t 
b vb 

^Y [ 2 d r Vc d u V vb (V ra Vsb (Kc V va - V ua V vc ) + V rb V sa (Ka Kc + Kc Ka)) + 



uvbc 



+ Q c ~ 2 5 r Kg 9„ Kb (Ka Kc + Kc Ka) (Ka Kc + Kc Ka) " 4 Kc Ka Kc Ka ■ (B8) 



S{r,(j) is a function of = 6 , and R a , which, being 3-coordinate-dependent, plays 
the role of an inertial potential for the gravitational field. 

Due to Eq.(3.14), the function <S(r, <r) appears in the Dirac Hamiltonian in the form 
f d 3 cr (1 + n) S (t,<t), so that we need the functional derivatives of this quantity. 
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a. The Functional Derivative with Respect to R a 



I 



dVi [1 + n(r, <7i) 



5S(t, S x ) 
5R a (T, a) 



+ 



= 2 (V E 

rsa 

<9 r ra (V m Ka [2 7aa _1 <9 S - (2 7 Sa 7^ - 5 d i) d s R b 

b 

+ Too Ka <9 S Ka + 7 5a Ka Kb Kb V Ka) ~ 

t)b 

- (1 + n) (V r aKa [2 7 aa(-0" 1 9 r 9 s + 30- 1 9 r 00- 1 9 s 0) + 

+ 2 0" 1 <9 r ^ (2 7 , a 7 - ba - 5 a5 ) <9 S i? 5 + 

6 

+ X] ( 2 7sa ~ ^ d r d s R~ b + Y^ t 2 75a (<*ac ~ Taa Tea) ~ Taa he] 9 r #b <9 S 
b be 

- 2(f)~ 1 d r (j) >Jaa Ka 8 S V ra + ^ 7ab Ka Kb Kb 9„ K 

vb 

+ ^ [(2 7aa 7Sa - 5 g5 ) <9 S (Ko Ko) + 2 ^ ( 76a Tab Taa Tbb) Ko Kb Kb d v K- 



+ 



d r R b - 



- 9 r (% a Ko <9s Ka + ^ 7afe ^ 9 « Vga ) + 

ub 

<9 r Kc <9n Kb (Ko Kb (Kc Ka - Ka Kc) + Kb Ka (Ka Kc + Kc Ka)) + 

uvbc 

+ (Taa - Tab + Tac) Qfe d r Kb <9 U Kb ((Ka Kc + Kc Ka) (Ka Kc + Kc Ka) ~ 

- 4KcKaKcKa)] )])(r,a). (B9) 
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b. The Functional Derivative with Respect to <ft = (j) 6 



I 



d 3 a x [1 + n(r, a x ) 



6S(t, &i) 
5 0(r, a) 



= 2 (0 (1 + n) [ ~ 8 (V 1 9r d s (V ra V sa ) + 

rsa 

+ V ra V sa {r l drd s <P-2r l d r <pY j r h adsRi)) + 

+ Vra V sa (2 ^ 7 - ba d r d s R b - ^ (2 75a 7c - a + <S fe ) d r R b d s ife) - 
b be 

- 2 E [V.« -9s Ka - 2 <9 S (Ka V,«)] + 

5 

+ Yl Tbb Vra Vvb Vsb dv Vsa ) 9r Rb ~ 9r 9s ( Vra Vsa ^ + 
vb 

+ \Y.[ 2dr Vsc du Vvb { Vra Vsb ( Vuc v ™ - Vua v ^ + Vrb Vsa ( Vua Vvc + Vuc v ™}) + 

uvbc 

+ Ql Qc 2 d r V sb d u V vb ((Vra V sc + V rc V sa ) (V ua V vc + V uc V va ) - 4 V rc V sa V uc V 

- Yl dr 71 Q° 2 [ 2 V ™ V ™ (4 0" 1 -9 S - ^ Tba 9s Rb) + O s (V ra V sa )] ) (t, B) . 

rsa b 

c. The Functional Derivative with Respect to 9 l 



(BIO) 



= -[0 2 ^Q" 2 9 r n(2 

rsa 

OV 



V ra V sa ,_i a 



- 2 



Y [ Tba Vra Ssv + Y Tbb Vva Vrb Vsb dv R ~ b + 

vb b 



dV M 



uvbc 
1 



+ Y ~~Q§T du Vvb Vra Vsb ( Vuc Vva ~ Vua Vvc ) + Vrb Vsa ( Vua Vvc + Vuc Vv < 



+ 



+ Qc 2 d u Vvb [(v ra v sc + V rc V sa ) (V ua v vc + V uc V va ) - 4 v rc v sa v uc V va ] ) 



r, 



uvbc 
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-2^ 1 9 r 9 s + 60~ 1 9 r 00~ 1 9 s + 



+ E ( 2 75a 7ca - <fe) <9 r & # 

d V ra V vb V sb 



be 

2 E ->«. 



90* 



Ka <9 r + 



s&6 



a. K ra + 



+ 2V ra ((f) 1 <9 r - E 7ca <9r #c) <9s #s - <9 S V ra d r Rij + 

c 

+ E 756 ^ rfe Vsb [ dr dv R ~ b + 

vb 

+ 2 9 r 0-E Tea d r Rc) & v + d r (Ka V rb V sb ) d v i? 5 ) 



+ 



+ E [^Kc^Kfe 

uvbc 



d 

dv 

06 



Ka V sb (Kc Ka - Ka Kc) + V rb V sa (V ua V vc + V uc Ka)) 
^(2 {r l dr<p-Y,r b adrR-b)d u V vb 



V ra V sb (V uc V va - V ua Kc) + V rb V sa (V ua V vc + V uc V va )\ + 
+ 8 r ]d u V vb (V ra V sb (Kc V va - V ua Kc) + V rb V sa (Ka Kc + Kc K 



uvbc 



- d r V sb d u V vb 



I ((Ka Kc + Kc Ka) (Ka Kc + Kc Ka) ~ 4 Kc Ka Kc Ka) ~ 
(2 (0- 1 d r - E (Tba - lib + 75c) d r R b ) d u V vb 

5 

((Ka Kc + Kc Ka) (Ka Kc + Kc Ka) ~ 4 Kc Ka Kc Ka) + 
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dVsi 
d6 i 



+ 



+ d r 



d u V vb ((Ka Kc + Kc Ka) (Ka Kc + Kc Ka) ~ 4 Kc Ka Kc Ka)] )] ) 



T,<7 



(Bll) 
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3. The Laplace-Beltrami operator A, the Function T and their Functional Deriva- 
tives 



The Laplace-Beltrami operator associated with the 3-metric 3 g rs (det g rs = 1) appearing 
in Eq.(3.43) is 

A = d r d 8 ) = 3 g rS 3 V r 3 V S = d r ( Qa 2 Va V sa d s ) = 

a 

= E Q* 2 V ra Vsa dr ds ~ ( 2 Vra Vsa X- %a dr Rl ~ dr {Vra Vsa) ) ds ] • (B12) 

a I 

From Eqs.(3.13) and (2.10) we have 
T(r,a) = ^9 r ( 3 e^V s3 9™R 3 fc-9 s 3 y)(r, ff > 

r suv 

= ~J2 d r (/ Qa 2 [ 2 V sa (4 0" 1 9 S - ]T Tba 9 S i2 6 ) + 9 S (Ka V*,)] ) (r, B) = 



= - [<^ 2 E 2 (2 K-a v; a [4 <9 r 8 S + 0" 1 «9 r 0" 1 d s 0) - 

rsa 

~ Y,^a{d r d s R l + 2(hr l d s <p-Y J l^d s R 5 )d r R-^\ + 

+ 2 3 r (Ka V sa ) (5 0- 1 9, - 2 ^ 7ba 9 S i2 6 ) + d r d s (V ra Vsa))] (r, <?) = 

6 

= 7 (-80A0 + T 1 )(t,5% 

rx(r, a) = - [0 2 5] g~ 2 (2 v; a v sa [a 0- 1 d r 0- 1 d s - 

rsa 

- E 75a (d r d s R- b + 2 (0 _1 <9 S - 7g a <9 S i2g) <9 r i2 6 ) + 

5 c 

+ 2 d r (V ra V sa ) (0- 1 d s - 2 ]T 75a d s Ri) + d r d s (V ra V sa ))] (r, a). 



(B13) 



In the Dirac Hamiltonian (3.48) there is the term — J d 3 a nT (r,a). To eval- 
uate its functional derivatives we use the form T = — 80A0 + T\. First we evaluate 



the functional derivatives of the term -^-q / d 3 a 



nd> A i 



16wG 



Jd 3 a nT (r,a) 



16ttG 



Jdh 



n( - 80A0 + T! 



r, er) and then of the quantity 
(r,a). 
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a. The Functional Derivative with Respect to R a 



I 



n d> A 



- cTT^ = 2 E ^ ^ 2 V ™ ^« ^ ^ ^ (" ^ d * * + 9 r n )] ( T > *)> 

a V T > <T J rsa 

(B14) 



/ dVi n(r, <n AB V - = " 8 / rf ai a p r ^ + 

y di2a(r,<r) y 5i? s (r,cr) 

+ 2 [0 2 ^ 7sa Q~ 2 V; a K a ( ( 9 r ^n + 29 r n0- 1 ^0 + 8n0- 1 9 r 00- 1 9 s 0)](r,(T) = 

rsa 

= 2 [0 2 ^ Taa Q- 2 V ra V sa [d r d s n-6 <T 1 d s <9 r n)] (r, a). (B15) 



( r , ^i) 



6. T/ie Functional Derivative with Respect to 4> 



d 3 ai 



5 nci> A oi> 



<5 0(r, <r) 

2nA0 + 2 Qa 2 V ra V sa d r nd s <j) + (f)An\ (r,a). 



(B16) 



.,/_ A ST ( T > gl) _ , / 

5 0(r, a) 



(i 3 <Ti n(r, 



-8 / dVi 



5 0(r, a) 



+ 



+ 2(8nA0 + 0g- 2 9 r n[9 s (\/ ra V; a ) + 

+ 2 \/ m V sa (4 0- 1 9, - £ 75a 5 S i2 6 )] ) (r, a) = 

= -2(0^g- 2 [4\/ m \/ sa( 9 r( 9 s n- 

rsa 

- 3 (2 V ra V sa ^ 75a 9 r R~ b - r (V ra Ka)) 9 S 71 ) (r, <j). 



(B17) 
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c. The Functional Derivative with Respect to 6 l 



( r , *i) 



E^_2 ^ Kg Kg 



n i 



~ l d r d> + d r n 



4 a 



(B18) 



/ 



<i 3 (Ti n(r, (Ti 



(JT(r,^) 



n (p A 



2 ^ Vsa 



d r d s n — 6d r n(j) 1 <9 S j(r, a). (B19) 



4. The Function 3 ,R. 



As shown in Ref.[6], for the evaluation of 3 R (a function of 9 l and R a ), appearing in Eqs. 
(3.43) and (3.44), we need the following results (see Eqs.(A25) of Ref.[2] and Eq.(2.10) for 
the spin connection 3 u; r ( a )) 



'R[6 n , 0, Ra] = J2 3e e W( b )( c ) 3(i U %) ^-( c ) = 3e 



abrs 



+ ^2 ( 3 ^r(a)(e) 3 ^s(e)(fe) - 3 ^s(a)(e) 3 ^r(e)(o)J 
e 

= (j)- 5 (-8A0 + 3 £0) =<P~ % (S + T), 



4 
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= Y,Qa 2 {-drd s {y ra v sa ) + 

rsa 

+ 2 E [r ba (2 & (K« Vsa) - v sa d s Ka) - ^ llb v ra v vb v sb d v v sa \ d r R- b - 

b vb 

- V ra V sa ^2 \d r Rbd s Ri + 2rba{^ lcad s Rc d r R- b - d r d s R- b y^ + 

b c 
[ 2d - V ™ ^ V * Vsb (Kc Ka - Ka Kc) + V rb Vsa (Kc Kc + Kc Ka)) + 



4 

uvbc 



+ Ql Qf 9r V sb O u V vb ((Ka V sc + V rc Vsa) (Ka Kc + Kc Ka) " 4 Kc Ka Kc Ka) ] ) ■ 

(B20) 



5. The Functional Derivatives of M. 



We need the derivatives with respect 9\ (ft and i? s of the mass density M.{r, a), appear- 
ing in the combination J d 3 a (1 + n) SA (t, a) in the Dirac Hamiltonian (3.48),given in 
Eqs.(3.37). They are 

^^(^(r))^ ((1 + n) 



Sj™ 2 C 2 + <ft~ 4 Ersa Qa 2 Vra Ka (/^M - f A ±r ) (^(r) - f A ±s ) 

+ (l + n(r, a)) (0- 1 / 3 

[ 1 n2 19 Ka Ka r s 1 n _ 2 n -2 g Ka Ka Kfc Kb p „ 

rsa abrsuv 

- ^ E ^ (2 - E ^ E ^ »* 5 m c(a, tf(r))) 

rsan m i 

5 sn J2 Qj Vj d n c(a, fjj{r))\ ) (r, a), 



(B21) 
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= -2^5 3 ( ( x,7/ i (r))r ?i ((l + n) 

i 

<t>~ 5 Ersa Qa 2 VraV sa (k ir {r) - f ^Lr) (W) - f 
^ C 2 + 0- 4 Qa 2 Ka Ka " f A ±r ) (^(r) ~ ^M ±s ) 

- (l + n(r,*)) (r 1/2 

[- ^ <5a V ra V sa 7T^ 7T* + — ^ Q~ 2 Q^ 2 Ka Ka Kfe Kb -^ru ^sj; - 
rsa abrsuv 

v ™ v » ( 2 < - E ^ rm E & »* 9 ™ c ^> ^( r ))) 

rsan m i 



(B22) 



E rsa 7saQ- 2 KaKa («ir(r) - f A ±r ) (/^(r) - f A ±s ) 

/ m 2 c2 + _ 4 E ^ g-2 y ra Vga ^. r(r) _^ lr ) (/^(t) - f A ±s ) 
+ (l + n(r, a)) 

7sa ^ ^™ ^ ^ _ 2c E ^ 7sa + 7sb ) 2 ^ ^ ra Vsa Vub Vvb Fru Fsv ~ 

rsa abrsuv 

- - c J2 ^ v ™ v «> ( 2 n ± - E E & ^ ^ ^( r ))) 

rsan m i 

3 

(B23) 
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APPENDIX C: THE QUANTITIES OF APPENDIX B IN THE 3- 
ORTHOGONAL GAUGES. 

In this Appendix we give the restriction of the quantities evaluated in Appendix B the 
the 3-orthogonal gauges. 

We shall use the notation introduced after Eq.(6.1). 

The functional derivatives of Eqs.(B2) - (B7) become 



n {a )\ei=o(r, a) = <p 2 (r 



tabi Qb 



Q»Q: 1 -Q.Q; 

tabi Q a 



^ , ~^ ^ 5a ~ 7a ) 9b R ~ c n i e) + 

b+a i [QbQa 1 - QaQb 1 ) 5 

+ Q- 1 (f d a H + ^2 ( 75a d a U- b - d a R- b U- b ) + M {a) )] (r, a) « 0. 



W( ) + 7l( (4 0- 1 9 a 0-^75a5 aJ R o ))](r,a). 



(CI) 



(C2) 



J n ( «)(r, *i) l*=o « 6 [0 3 £ fi (o) Q" 1 4 ttJ (r, a). (C3) 



/ 



d (Ti "(o)(t, o-i ) I ^, — I" If- o 



^iab Qb 



5 nf" 1 (r, a) 

.1 -r 



a^6 



QbQ^ -QaQb' 



-<f> 1 d b (f) + ^2rbadb R b ) n {a ) - d b n {a ) ) (r, a). 



(C4) 
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SO^r, a) 
-(f)' 2 (t, a) E d r n {a) (T,a) 

ra 

[Qa 1 E ( 75 « V ^ ra + Tbr V ^ ar ) U ~ b + 



1^=0 — 



+ e(ee 



c -abtQ 



-1 



t b^a 



Qb Qa 1 ~ Qa Q\ 



+ ZZ 



e 6cj (5 rc V(i) ab — 8 ra V(j) c j,) 



QcQb 1 - QbQc 1 



TV 



(0) 



(r, <?) + 



+ 2 ( r >^) 5^ n («)( r '^) 

a 

1 ( E v w™ V" d * n + E (> ( 2 ^ ^ + E ^ ^ n s 

r fe c 

- d^IIg) +E TbadrRb + Mr] ~ 

b 

- E V W<* Tbb [9b n 6 - (2 0" 1 <9 fe + E 7c-a «9 b i? c -) n- b ] ) - 

bb c 

-EE ^ c ^ - ^ *W) * + 



rj b c ^ fe ^Qb 1 - QbQc 



+ E (EE 
+ EE 



Qfc Qa 1 - QaQ 
e&ci Qa 1 QbQc 1 



jr t b^a ^"^a -su-^b 



-(2 1 8 r + E 7ca <9r #c) V(i) r6 + <U (V(i)jt + B {i)jt )\ + 



c^b 



QcQ b —QbQi 



J (Src V(i)ab - &ra V(i)cb) [2 1 <9 r + 



QcQ^ + QbQ' 1 



^ ^ ( Iba Ibb Ibc {ibb IbcJ , j , , s_ 



<9 r i? b 
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1-1/3 



E ( E ^ n ( fl ) ®° 1 E ^ + 7 ^ V ^ U ~ b 



^{^Qb 1 V (i)rb ° (a)(6) + ^ QbQr- 1 V (i)ab ° ( 

bj^a bj^r 

tart Qr 1 (Y(i)jt + B(i)j t ) tjcd Qc Q d 1 



r)(b) 



EE 



tj c^d 
^3 



0-(c)(d)) 



+ 



8ttG 



1 Qa Qr 1 

<9a n (o ) Qa 1 ^2 QbQc 1 V (i)cb ^(6)(c)) ] (t, 5) 



+ 



+ 



b^c 



+ 



+ [r 1/3 E "w (q- 1 E [^(06- - E ^ ^ n 5) 

a 6 6 

+ E ^ v w a ~ Tbb v ^ dh u ~ b + 

+ E ( Tba v ^ ba + Tbb v ^ ab "> ^ + E 7ea ^ ^) n & + ^(o&o ^ 



+ 



+ 



8ttG 







/ 2 

& (6)(c) - V(i)c6 <9a CT (b)(c) + ~ ( V(i)a6 1 <9 C 



- V( i)c6 1 <9 a 0) <T(6)( C ) - E i^ba ~ Tbb + lie) (V(i)ab C R~ h - V {i)cb 8 a R~ h ) <7( 6 )( c )) + 

b 

+ y y eabt Q ' bl {Vm I Bm) ej T Qc — ( - r 1 9 b + y Tba d b R- b ) *< em 



tj b=£a,c=£d 



/ #<71 E 

^(a 



_ <JW(o)(r,5i) 



5n s fr, a) 



|6» 1 =0 



_<r 2 E 1 ("w 

a 

7aa<9 a n (a) j] (r, a). 



7„ a ^2 1 (9 a + E 75a <9a #b) - <9 a 
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J dji y M(T|Jl) gifefoff) '^° = 
= <j>~ 2 (r,a) d r n {a) (r : a) Q' 1 5 ra Tl a - 

ra 

- 2 £ ( 7a „ - 7at ) £ ^g^' -y (r, a") + 

a 

q- 1 (a a n a - (2 0- 1 & + 75a 4 i? 5 ) n a - 

5 

- Taa [<P 6 da^ + Y, (iba 9 a ^ - d a R~ b U~ b ) + M a ] ) - 

- Taa J2 E 



b^a i 



+ ?££ 

i b^a 

<r<»] (r, 3) 



^abi Q a 
(QbQa 1 ~ QaQb 1 ) 2 



2 (7aa - Tab) 1 <9 b + ^ (7 aa 7bb ~ Tab Tba) & R~ b 



E 1 "(«) n " + 4^ ^ ^aa ~ 7ab 

b^a 



<9b^(a) C(o)(6) 



+ 



-1/3 



QbQa'-QaQ, 

«(«) (Qa 1 [da ^ - {\ 0^ d a + £ 7 ba 5a i2 5 ) II a - 

a ft 

- 7aa (09 a ^ + J](7 5a 9 a n 5 -a a i? 5 n- 6 ) + A4 a )] + 

6 

+ 8^G^ ^ ^ pa9 6 (7 (a)(b) + 



+ 



b^a 
2 



+ ^(7aa g (7aa 7ab) 



QbQa'-QaQ, 



-) 0-^,0 



- E 



7aa (7ba ~ 7bb) [gb Q a 1 + Qa Qfc j + 2 (7aa 7bb ~ 7ab 7ba) Qb Qg 1 

QbQa 1 ~ Qa Qb 1 



8 b Rlj <T(a)(b) 



The function S(t, a) of Eqs.(B8) becomes 
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S\ d i =0 = <p 2 ^ Qa 2 (Y [ Yl ( 2 Tha 75a ~ ^) da R ~ b da R ~° 

a fee 

- 4 Tba 0" 1 d a 8 a R- b ] + 8 d a 0) 2 ) , 
and its functional derivatives (B9), (BIO) and (Bll) become 



/ 



dVi [1 + n(r, CTi)] 



£«S(r, gi) 
SR a (r, a) 



le«=o 



= 2 (0 2 ^ g^ 2 [9 a n (2 7 «a 0" 1 9 a - Y ( 2 7aa 75a ~ d a R~ b ) 

a b 

- (1 + n) (2 Taa ( - 0" 1 d 2 a + 3 9 a 0) 2 ) + 

+ ]T (2 7 aa 75a " hi) i% + 2 0" 1 9 a d a R~ b ) + 
5 

+ Y i^Tba{^ac - laalca) ~ laak^j d a R- b d a R^j )(t,<t). 



be 



(C8) 



(C9) 



d 3 <7i [1 +n(r, ai)] 



SS(t, <7\) 



=0 — 



(5 0(r, <r) 



+ (1 + n) (8(0- 1 9 2 0-20" 1 9 a 0^75a5a J R o )-2^75a5 2 i? E + 

5 5 

+ ]T (2 75a 7c"a + 55c-) ^a i?5 ^a i?c") 1 ) (t, 5) . 



be 



rfVi [1 +n(r, cti)] 



J¥(r, a) 



(CIO) 



= -2 [0 2 Q a 2 V(i) ra (d a n((j) 1 <9 r - Y Tba d r R- b ) + d r n((j) 1 d a - 7 5r <9 a i? 5 ) - 

5 5 

- (1 + n) [ 2 (3 0- 1 ,9 a 00- 1 ( 9 r 0-0- 1 ( 9 a( 9 r 0) + ^(75a + 75r)«9a«9 rJ R- fe + 

5 

+ 2 Y ( 75a 5a ( t )d r R b + Tbr 0" 1 <9 a R~ b ) ~ 

5 

- ^(2 7 5r 7 c - a + 5 5c -)«9 aJ R 5( 9 rJ R c -])](r, ( 7). (Cll) 
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The Laplace-Beltrami operator of Eq.(B12) becomes 



A|, i=0 = Qa 2 k 2 - 2 r ba d a R- b d a 



(C12) 



The functions T(r, a) and Ti(t, a) of Eq.(B13) become 



T(r, a) 



-2 [0 2 £ g- 2 ( 4 9 2 + 9 a 0) 2 " - 

a 

- E Tba [dl Ri + 2 (5 0" 1 9 a - 7c"a d a i? c -) 6> a ) ] (r, a) , 

b 5 

71(^)1^=0 = -2 [0 2 £ g- 2 (4(0- 1 9 a 0) 2 - 

a 

- E Tte Ri + 2 (<T 1 5 a - E 7c-a 9 a R 5 ) d a R- b ] )] (r, a). (C13) 



The functional derivatives (B14)-(B19) become 



/ 



nd> A< 



=0 — 



2 2 ( E G« 2 7™ ^ d * <t> (n 0" 1 9 a + d a n)) (r, a), 



(C14) 



rf 3 a"i n(r, c?i) 



5R a (T, a) 



= 2 [0 2 E 7aa g~ 2 (d a 2 n - 6 <T 1 «9 a d a n)] (r, a) 



(C15) 



/ 



n0 A 



^ 3 o"i 7—7 — =^ 1^=0 - 

0(r, cr) 

2n A| e , =o + 2 E ga 2 9 a n9 a + 0A| ei=o 



n 



(C16) 



/ 



d <7i n(r, <?i) 



^^T(r,5i). 



5 0(r, a) 



|0»=O 



-4 [0 E g- 2 (2 9 2 n - 3 E 75a 5 a i% d a n)] (r, a) 



(C17) 
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n <j) A . 
J¥(r, a) 



(T, ^1 ) 



+ drucj)' 1 d a (j) + d a n(j)~ l d r (j) )(r, a), 



(C18) 



/ 



<i 3 <Ti n(r, cti) 



<JT(r,5i) 



^"1^=0 



= -2 J0 2 Q a 2 V(i )ra (d r d a n - 3 (<9 r n0 1 d a (j) + d a n(f) 1 <9 r 0)) (r,a). 

ra 

The scalar 3-curvature of Eq.(B20) becomes 

3 %,a)|, i=0 = £ (g a - 2 ^ [2 75a 9 a 2 J R 5 -(9 aJ R 5 ) 2 - 
- 2 7 5a 7 5a 9 a i?g <9 a i? 5 ) (r, a) . 



(C19) 



(C20) 



Finally the functional derivatives (B21), (B22) and (B23) of the mass density become 



da, (l + nfo^)) |^o = 

= \Ys^Mr))m{{l + n) 

i 

0" 4 Ersa Qa 2 (V(i)ra + ^ra V( i)sa ) (/^(t) - & A ±r ) (^ s (t) - & A A 

^jmj C* + £ a g- 2 (^ a (r) - f A ±a ) 2 

+ (l + n(r, s )) (0" 1/3 ~ c [ Ql V{i)ra Ssa7C r ± < + ^ ' ^ V(i)r« Kb 

ars afer 

- (V(i)ra + V (l)sa <5 m ) (2 TT^ - ^ ^ £ Q< ^ 9 m C («7, 7fc(r))) 

arsn m i 



)(r,<?) 



+ 



(C21) 
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= -2^5 3 (a,77,(r))77,((l + n) 

i 

^mf + 0~* Ea q-2 ^ m(t) -%A ±a ) 

- (l + n(r, a)) (> 1/2 [± £ Q 2 5 m 5 sa ^ < + 1 ]T Q~ 2 Q^ 2 F afe F afe 

ars afe 

- ^ E Ql 5 ™ Ssa (l ^ - J] 5™ £ 7ft 9 m C (<7, T^r))) 

arsn m i 

<5 sn ^ #j Vj dn c{a, ^(r))] ) (r, a), 



(C22) 



-X)5 3 (a,^(r))i7i ((1 + n) 



=)(r,<?) 



^/m 2 C 2 + 0-4 Ea g- 2 ^ w(r) _ fit A±a y 

+ (l + n(r, <?)) (0- 1 / 3 [1 £ 7sa Ql 5 ra 5 sa n r ± vri - 1 ^ ( 7sa + 7sb ) Q~ 2 Q fe ~ 2 F a6 F afe - 

ars " afe 

- \ 7a« <*ra 5 S a (2 7^ - £ 5™ £ Qj Vi d m C (<7, 7fc(r))) 

arsn m j 

j 
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